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Summary 

A four-cylinder-Stirling-engine, transient-engine- 
simulation computer program is presented. The program 
is intended for controls analysis. The associated engine 
model has been simplified to shorten computer 
calculation time. The model includes engine mechanical 
drive dynamics and vehicle load effects. The computer 
program also includes subroutines that allow (1) 
acceleration of the engine by addition of hydrogen to the 
system and (2) braking of the engine by short circuiting of 
the working spaces. Subroutines to calculate degraded 
engine performance (e.g., due to piston ring and piston 
rod leakage) are provided. Input data required to run the 
program are described and flow charts are provided. The 
program is modular to allow easy modification of 
individual routines. Examples of steady-state and 
transient results are presented. 


Introduction 

A Stirling engine is a mechanical device that operates 
on a closed regenerative cycle. Cycle expansion and 
compression of working fluid is done at different 
temperature levels. Flow is controlled by volume 
changes. Heat is supplied continuously and externally 
and extracted with a net conversion of heat to work. In its 
simplest form a Stirling process consists of two pistons, 
one operating in a hot space and one in a cold space. A 
working fluid between the pistons moves continuously 
back and forth within a working space between the hot 
and cold spaces, where it is continuously heated and 
cooled. The fluid passes through a regenerator, which is a 
heat storage device. The regenerator stores or gives off 
heat depending on the temperature of the fluid passing 
through it. The result of the process is the Stirling cycle. 
Stirling engines are used for stationary power generation, 
heat pumps, refrigeration systems, and artificial hearts. 
Another application is the automotive Stirling engine. In 
this application there are four pistons, all interconnected 
by four working spaces (fig. 1). This report presents a 
computer simulation of an automotive Stirling engine 
that can be used for controls research. 

Stirling engine modeling is usually done for the 
purpose of predicting engine performance. The models 
are, of necessity, complex since Stirling engine 
performance (e.g., efficiency and fuel consumption) is 
quite sensitive to factors such as pressure drop within the 
engine. In general a Stirling engine performance model 
includes two pistons and the working space between 
them. The pistons are positioned as a function of 
time and are 90° out of phase. Control volumes are 
used to segment the working space into its 
components - expansion space, heater, regenerator, 


cooler, and compression space. More than one control 
volume can be used to describe a component. Various 
heat transfer paths are usually included. Although many 
Stirling engines have more than one working space, the 
models used to predict performance generally do so for a 
single working space, and the resultant power generated 
is multiplied by the number of working spaces. One such 
model has been developed by Tew, Jeffries, and Miao 
and is described in reference 1. It contains 13 control 
volumes in the working space. The model calculates flow 
resistances and heat transfer coefficients as well as 
performance over a single engine cycle. Daniele and 
Lorenzo (ref. 2) have reduced the number of control 
volumes in the working space to seven, one each for the 
expansion space, heater, cooler, and compression space 
and three for the regenerator. Also, a metal volume has 
been associated with each regenerator gas volume. 
Average flow resistances and heat transfer coefficients 
are derived from the Tew model. Within each gas volume 
the continuity and energy equations are integrated and a 
simplified, first-order momentum term (flow resistance) 
is calculated. Upwind differencing (ref. 3) is used to 
calculate the interface volume temperatures for use in the 
energy equation. The resultant model has 17 state 
variables and uses a backward-difference integration 
scheme for problem solution. Results generated by that 
model compare well with experimental power and torque 
data over the entire speed and pressure ranges of the 
engine. 

The aforementioned modeling approach produces 
acceptable steady-state (power and torque) predictions 
when working spaces in the engine are isolated. However, 
proposed control schemes for four-cylinder Stirling 
engines require that the working spaces not remain 
isolated. One such scheme is engine braking by short 
circuiting. In this case, all four working spaces are 
connected to high- and low-pressure manifolds. Flow 
enters or leaves a working space depending on the 
pressure differential and the orientation of the check 
valves between the pressure manifolds and the working 
spaces. A complete four-cylinder model is needed to 
analyze this controls problem. 

One such four-cylinder controls model is described in 
reference 2. In that model, each working space model is 
further simplified to reduce computer run time. This 
simplification involves reducing the number of control 
volumes in a working space from seven to three - one for 
the expansion space and the heater, one for the 
regenerator, and one for the cooler and the compression 
space. It is also assumed (1) that the temperature within a 
control volume remains constant and (2) that the gas and 
regenerator mesh temperatures are equal. This reduces 
the number of state variables in a working space to three 
(or 12 for the whole engine). Results from this simplified 
model are presented in references 2 and 4 and agree well 



with steady-state experimental data. The model has been 
used to predict the steady-state performance of the engine 
with short circuiting and with piston ring leakage. 
Reference 4 describes the use of the simulation to study 
supply transients with and without piston rod leakage. 

This report documents the simplified four-cylinder 
model and its implementation. A users manual is given 
on how to run the simulation. Both steady-state and 
transient results are presented and a printout is provided 
for a test case. Finally flow charts are given for the 
subroutines and the overall simulation. 


Model Description 

A schematic of the four-cylinder-Stirling-engine 
controls model is shown in figure 1 . The four cylinders 
are shown interconnected by the four working spaces. 
Each working space contains three volumes - one for the 
expansion space and the heater, one for the regenerator, 
and one for the cooler and the compression space. Each 
volume is assumed to be at constant temperature. 
Variations in volume temperature can be scheduled as 
a function of engine speed. Variable names shown in fig- 
ure 1 correspond to the computer coding. A complete 
symbols list is given in appendix A. 

Besides the engine thermodynamic elements, three 
other system models are shown in figure 1 - a drive 
geometry model, a working-fluid supply model, and a 
short-circuiting model. Torque is calculated as a function 
of differential forces on the pistons and summed in the 
drive geometry model. Losses due to auxiliaries, load 
effects, etc., are included. The net torque is then 
integrated once to give engine speed and again to give 
crank angle. The remaining models shown in figure 1 are 
associated with control schemes for accelerating and 
decelerating the engine. The acceleration system is 
indicated by alternate long- and short -dashed lines and is 
described in appendix B, Engine deceleration is 
accomplished by short circuiting the working spaces. This 
deceleration system is indicated by short -dashed lines. 
This system is also described in appendix B. 

A schematic of the drive dynamics is shown in figure 2. 
Differential forces on the pistons are translated into 
torque through the vehicle drive geometry and summed 
to form total torque (TORQT). Torque due to engine 
friction is subtracted to form brake torque. This is 
available to drive the auxiliaries and the vehicle load. The 
vehicle inertia and the gear ratio are used to compute the 
effective load. The summation of torques is integrated to 
give engine speed (PSIDT) and again to give crank angle 
(PSI). The crank angle is used to generate piston position 
(by using the crank geometry). The model can be run in 
this manner, or piston position can be input as a function 


of time with torque and cycle performance calculated at 
constant speed. 

Relative to the performance models, the controls 
model is simple, yet it contains the essential elements for 
controls analysis. It has been shown in reference 2 that 
simulation predictions compare well with experimental 
data. However, it should be noted that heat flow 
dynamics are not included in the model. These effects 
may, in some instances, be significant and thus necessi- 
tate expansion of the model. 


Users Manual 

Simulation Flow Diagram 

The overall simulation structure is shown in figure 3. 
Run conditions are set in the main program (MAINSE). 
Subroutine ICSTUP is then called to calculate initial 
conditions. SHORT2 and PISTN3 are called from 
ICSTUP to calculate the initial conditions for short 
circuiting and initial piston positions, respectively. 
ICSTUP then calls FOURW2, which is the integration 
subroutine. FOURW2 handles the incrementing of time, 
data output, and run termination. FOURW2 calls the 
Stirling engine simulation subroutine FWS3V2 to 
generate the information it needs for the Jacobian matrix 
associated with the backward difference integration. 
FWS3V2 calls LOSSES to calculate engine friction and 
auxiliary losses and vehicle load effects; PISTN3, to 
calculate piston positions from crank angle and crank 
geometry; PLEAK2, to calculate leakage flow if piston 
ring leakage effects are considered; SHORT2, to 
calculate short-circuiting flows if engine braking is 
desired; and SUPPL2, to calculate the flow into the 
engine and the phasing of the injected flows. SUPPL2 
calls SUPLRK if rod leakage is considered. FWS3V2 then 
returns to FOURW2. 

In figure 3, the body of FOURW2 is shown within the 
dashed line. The subroutines enclosed are actually 
subroutines to FOURW2. Although they are all called by 
and return to FOURW2, they are shown as calling one 
another to illustrate the looping that takes place within 
the subroutine. Once a Jacobian matrix is calculated, 
FOURW2 calls DMINV for a matrix inversion or 
BROYFl for updates to a previously generated matrix. 
Subroutine BROYFl contains the Broyden update 
algorithm. This algorithm is used to try to eliminate the 
need to calculate new Jacobian matrices as the simulation 
moves away from an operating point. The algorithm tries 
to accomplish this by continually updating the inverted 
Jacobian matrix. The reason for using the algorithm is to 
try to shorten computer calculation time since generating 
Jacobian matrices and taking inverses is very time 
consuming. The use of this algorithm is discussed in more 
detail in appendix C. Convergence is then checked. If the 
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Figure 2. - Stirling engine drive dynamics. 



Figure 3. - Overall simulation structure. 


simulation is not converged, more passes are made 
through FOURW2 to generate better iteration guesses or 
a new Jacobian matrix. If converged, FOURW2 calls 
TRAP to do a trapezoidal integration on the product of 
pressure and volume in the expansion and compression 
spaces. TRAP returns to FOURW2, which then calls 
OUTPUT if a printout is desired at the current time. 
OUTPUT returns to FOURW2; then GUESF2 is called 
to predict the next set of state variables for the next time 
step. This is done until the desired maximum number of 
time steps has been run, at which time FOURW2 
terminates the run. Flow charts for all subroutines are 
given in appendix C. 

Program Setup 

The program setup is done in the main program 
MAINSE: Switches are set to indicate the type of 
transient desired, and engine geometry and flow data are 
input. Tables I to X indicate the required input as well as 
the options available. In table I engine geometric data are 
specified. The expansion-space dead volume (VOE) 
includes the heater volume; the compression-space dead 
volume (VOC) includes the cooler volume. The total 
regenerator volume excluding the volume of the mesh is 
denoted by VR. Table II lists the required heater and 
cooler wall temperatures. Table III lists the flow 
resistances between the volumes. All the constants are 
listed in table IV. Data for the implicit integration 
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method are listed in table V. Most of the integration 
settings should be kept the same. However, increasing 
MPAS and/or TOLPCQ sometimes is helpful for 
difficult convergence problems (not observed in this 
application). The simulation will continue to run if 
MPAS is exceeded but will output an error message 
indicating it has been exceeded. 

The different options available can be selected by a set 
of switches defined in table VI. ISS must be set equal to 
zero and MATRIX, equal to 1 . The rest of the switches 
are user definable. Table VII lists all the run conditions 
that must be specified. Table VIII defines the load 
characteristics. Table IX describes the cycle data -the 
number of cycles to be run and the number of 
integrations per cycle. Table X lists the required transient 
data. If no supply or short-circuiting transients are 
desired, the cycle start-and-stop conditions for these 
options should be set higher than the number of cycles 
desired. Table XI lists all the common blocks and the 
subroutines that contain them. This table is provided so 
that the user may more easily define data transfer 
between subroutines. 

Output Options 

The user may select from a number of output and 
debug options. If ICALC = 1, a short printout is 
specified. The output, where all symbols are defined in 
appendix A, is given in table XII. Outputs include heater 
temperature (TWH), cooler temperature (TWC), cycle 
mean pressure (CYCLPR), engine drive mode 
(NONENG), number of cycles to be run (NUMBCY), 
piston ring leakage area scalar (ALEAK), piston rod 
leakage area (ARLEAK), supply pressure (PSOURC), 
the time step to start the supply (ISUPST), and the time 
step to start short circuiting (ISHTST). No output values 
are given in this example since this output option was 
used for all the test cases. 

Next, results from the output data are printed out at 
the desired time steps. Power and torque calculations 
along with crank angle, mean pressure, and engine speed 
are output. ITRAN indicates the number of time steps 
taken up to the printout; and KWORK indicates the 
number of time steps in the calculation of the pressure- 
volume area. (Note that if KWORK = 101 , 100 time steps 
were taken, since the first time step printed out is the 
initial condition.) 

If ICALC = 0, a long printout is specified. An example 
of a long printout is shown in table XIII. At each time 
point a complete listing of all engine variables is given. 
Note that the output is arranged in rows of four (except 
for the overall parameters). The first row is for the first 
working space; the second, for the second working space; 
etc. The output includes pressure, temperature and 
volume data, flow data, convergence data (guess 
variables (VS) and errors (E)), the short-circuiting 


system, and finally, power and torque data. Data are 
given for this printout as an example; subsequent 
examples will make use of the short printout 
(ICALC = I). 

There is also a debug option to aid in solving problems 
that may occur in the simulation. This option is specified 
in MAINSE by NOBUG. If NOBUG = 1, the option is 
not used unless (1) there is a problem with convergence 
(the maximum number of allowable iteration passes 
(MPAS) is exceeded without convergence) or (2) there is a 
problem in generating a partial derivative for the 
Jacobian matrix. In both cases the debug option comes 
from FWS3V2. No printout of variable names is given. 
However, an indicator shows where and what subroutine 
the debug output comes from. For example, if “DEBUG 
PRINTOUT FROM FWS3V2 NUMBER 1“ is printed, 
the user can go to FWS3V2 and look for the comment 
that indicates debug printout number 1. Using the data 
printed and the FORTRAN listing of the variable names, 
the user should be able to debug his program. 

If the debug option is called (NOBUG = 0), a dump of 
all variables as they are calculated comes from FOURW2 
and FWS3V2. This option should be used judiciously 
because of the amount of data printed. Again, an 
indicator points to the source of the debug output to help 
the user find the source of a problem. A sample debug is 
shown in table XIV. The first part of the output comes 
from ICSTUP and gives all the initial conditions 
calculated. Next comes the heading “DEBUG OUTPUT 
FROM FOURW2 NUMBER 1.’’ Here the variable 
TIME is listed. Next “DEBUG PRINTOUT FROM 
FWS3V2 NUMBER 1” is printed. Here the variables are 
not labeled, so the user must trace through the 
FORTRAN for FWS3V2 to the proper area of code to 
begin debugging the problem. In the sample the second 
debug output from FWS3V2 is also given. The output 
“CONVERGENCE IN ERR VECTOR” is printed if 
NOBUG = 0 and convergence occurs. 


Output - Test Cases 

Two output test cases -a supply transient and a short- 
circuiting transient are given. Input and output data are 
provided for both cases. 

Supply Transient (100 Points per Cycle) 

A supply transient is used as the first test case. The 
assumed initial conditions for the engine are a mean 
pressure of 5 MPa and a speed of 2000 rpm. After five 
cycles, hydrogen is injected to accelerate the engine for 
1000 engine cycles. The supply pressure is 10 MPa and is 
constant. Leakage through piston rings and piston rods is 
not considered. 
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The FORTRAN input in MAINSE for this case is 
shown in table XV. The transient is defined by 
CYCLPR = 5.0, SPDRPM =2000.0, NCYSUP =5, 
NUMBCY = 1000, PSOURC=10.0, ALEAK = 0.0, and 
ARLEAK = 0.0. For this case the Broyden update 
algorithm is not used (IBRYTH = 0); the step size is set by 
the desired 100 integration points per cycle 
(NPTPCY = 100); the desired printout is at every 100 
points (IPRTOP = 100); the short printout option is 
desired (ICALC=1); and no short circuiting is used. 
(NCYSHT is set above NUMBCY, or at 50 000).) 

A sample output for the test case is shown in table 
XVI. Note that TIME is incremented in such a way that 
the specified number of integration points per cycle is 
obtained. Therefore the time step is variable as a function 
of engine speed. For the example presented, only parts of 
the printout are shown because of the large number of 
engine cycles simulated. Before a specified limit on CPU 
time was exceeded, 803 engine cycles (20.77 sec of engine 
time) were simulated. The simulation was run on an IBM 
370/3033 computer and used 20 minutes of CPU time. 
Thus the simulation takes about 1.5 seconds of CPU time 
per engine cycle. Figure 4 shows a plot of the results. 
Note that the engine mean pressure (PAVEMP) and 
output power (POWERT) rise quickly. Engine speed rises 
more slowly because of the effective inertia of the load. 

Supply Transient (25 Points per Cycle) 

The second test case is the same as the first case but has 
NPTPCY = 25 instead of 100. Sample output is presented 



0 5 10 15 20 

Time, sec 

Figure 4. - Supply transient (100 integration points per cycle). 


in table XVII. Note that KWORK = 26 since there are 25 
points per cycle. POSDEG is somewhat off after each 
cycle (270.2°) since the integration is coarser than for the 
100-points-per-cycle case. Also, at some time points the 
simulation does not converge. Since MPAS is set equal to 
20 in MAINSE, the simulation prints a debug output 
after 21 iteration passes. The message “ITERATION 
FAILURE 13 21” “DEBUG OUTPUT FROM FWS3V2, 
FLOW CONTINUES” contains the number of 
converged errors, 13 (this will vary), and MPAS plus 1, 
21. Note in the first debug printout that there are 13 
converged errors but 14 states. Error 7 is -.115E-3, 
which is greater than the specified error tolerance 
(TOLSS in MAINSE), — O.lE-3. This convergence 
failure occurs because this is a rather coarse integration 
(large step size). The simulation continues, however, and 
eventually recovers. 

Since the simulation fails to converge, there may be 
some question about the validity of the output. Figure 5 
shows a comparison of the 25- and 100-points-per-cycle 
cases. There is excellent agreement between the cases 
because the errors that do not converge are very close to 
the tolerance band. For the 25-points-per-cycle case, all 
1000 cycles were simulated in 14 minutes of CPU time on 
the IBM 370/3033 computer, or about 0.81 second of 
CPU time per cycle. This amounts to 25.33 seconds of 
engine run time. 

As a matter of interest, the 25-points-per-cycle case 
was rerun by using the Broyden update algorithm 



Time, sec 


Figure 5. - Comparison of supply transients for different integration 
step sizes. 
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(IBRYTH = 1). The case took about the same amount of 
CPU time for 1000 cycles as it did without the algorithm, 
the only difference being that two more iteration failures 
occurred. It may be that the number of calculations 
required for the algorithm to update the inverted matrix 
approaches that required for generating a new matrix and 
its inverse for low-order systems. The algorithm is 
supplied in this simulation package to allow for 
expansion of the simulation to include more volumes, the 
energy equation dynamics, and a control system. 

The increase in the time step (decrease in the number of 
integration points per cycle) is desirable to decrease the 
computer run time. However, as the integration becomes 
coarser, some iteration failures can occur. Good 
agreement may still result, depending on how close the 
unconverged error was to the tolerance band at the 
convergence failure. This is problem dependent. The 
same supply case was run with 50 integration points per 
cycle. No iteration failures occurred, and 1000 engine 
cycles took 18 minutes of 370/3033 CPU time, or 1.1 
seconds of CPU time per cycle. Again good agreement 
was obtained with the 100-points-per-cycle case. 


allows for user selection of the desired transient and for 
detailed study of all or part of the transient. Results can 
be obtained over many engine cycles (as in the test cases) 
or on an individual cycle basis (as was done in ref, 4). 


Concluding Remarks 

A four-cylinder-Stirling-engine computer program is 
presented. The program is intended for controls analysis. 
The associated engine model has been derived from a 
more complex four-cylinder engine model. The simpler 
controls model results in a decrease in computer 
calculation times. The model includes drive dynamics and 
vehicle load effects. The computer program also includes 
subroutines that allow simulation of control strategies 
such as acceleration of the engine by adding to the 
inventory of hydrogen in the system and braking of the 
engine by short circuiting between the working spaces. 
Subroutines are provided to calculate degraded engine 
performance due to piston ring leakage. 


Short-Circuit Transient 

The final test case is a short-circuit transient. This is set 
up in MAINSE by setting the short-circuit valve area 
(ASHTT), the cycle at which to start the short circuiting 
(NCYSHT), and the cycle at which to stop the short 
circuiting (NCYSHS) and by setting the supply start 
(NCYSUP) and the supply stop (NCYSTP) cycles greater 
than the maximum number of cycles (NUMBCY). The 
plenum pressures PSHTMX and PSHTMN are set 
internally in the program to 1 MPa above and 1 MPa 
below the mean pressure, respectively. The MAINSE 
statements for this case are given in table XVIII. The 
output for this case is shown in figure 6. Maximum and 
minimum plenum pressures are shown at the top of the 
figure. Their initial values are 16 and 14 MPa, 
respectively. At five cycles (126 time steps) a short- 
circuiting sequence is begun. The valve area (ASHTT) is 
set equal to zero initially so that the plenums come to 
maximum and minimum pressures of 18.9 and 11.7 MPa, 
respectively. Then, 500 time steps later (internally set in 
the program), the short-circuiting valve is opened. The 
plenum pressures respond immediately, as indicated. 
Total engine power drops from 56 to 36 MPa when the 
valve opens. Speed drops off more slowly because of 
inertia effects. The torque, after 18 seconds of engine run 
time, is still slightly negative. Hence, speed continues to 
decrease. Mean pressure in the engine stays fairly 
constant. 

The test cases show that the simulation is capable of 
producing transients that are representative of Stirling 
engine controls analysis. The input routine (MAINSE) 
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All input data required to run the program are 
described, and complete FORTRAN listings of the main 
program and all the subroutines are provided. Flow 
charts of the overall simulation and the subroutines are 
also given. The simulation is modular to allow for easy 
modification of the model. 

An input routine is provided in which the user specifies 
the transient to be run by setting switches and supplying 
the proper geometric and engine data. The output is also 
selected by the user by setting a switch in the input 
routine. A very detailed printout at the desired printout 
increment or a less detailed printout can be selected. The 
integration step size and the printout interval do not have 
to be the same. 


The simulation has been run to generate both steady- 
state and transient results. Power and torque predictions 
from the simulation compare well with experimental data 
over a wide range of engine speeds and pressures. The 
simulation is capable of predicting steady-state 
performance with piston ring leakage and short circuiting 
between the working spaces. Engine acceleration and 
deceleration data are presented, along with associated 
printouts that can be used as test cases. 

Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, February 26, 1982 
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Appendix A 
Symbols 


A 

AAN 

AD 

ADJCMP 

ADJEXP 

ADM 
ADR 
ADRM 
ADSl 
ADS2 I 
ADS3 f 
ADS4 J 
AGDNl 1 
AGDN2 ^ 
AGDN3 J 

AIE 

AINERT 

AIVEH 

AIWHEL 

AKGl I 

AKG2 V 

AKG3 J 

ALEAK 

ALPHA 

AMPLIT 

AO 

APSI 

APTS 

AR 

AREA 

AREAO 

ARLEAK 

AS 

ASHTT 

RIGA 

BIGDEL 

BIGNUM 

CCCHK 


piston ring leakage area, cm2 
problem order 
piston area, cm2 

adjustment to compression-space 
temperature due to speed, K 
adjustment to expansion space 
temperature due to speed, K 
piston area, m2 
piston rod area, cm2 
piston rod area, m2 

slot area in piston rods 1 to 4 for 
working-fluid supply, cm2 

adjustments to cycle pressures 2, 3, and 
4, kg/(m2 N) 

engine moment of inertia, N-m/sec2 
total effective moment of inertia, 
N-m/sec2 

equivalent load moment of inertia, 
N-m/sec2 

wheel moment of inertia, N-m/sec2 

adjustment to cycle pressure, N-m/kg 

piston ring leakage area scalar 
crank angle lag, deg 
half of piston stroke, cm 
piston ring leakage area, cm2 
number of crank rotations 
number of integrations per cycle 
piston rod area, cm2 
trapezoidal integration area, N-m 
summed trapezoidal integration 
area, N-m 

piston rod leakage area, cm2 
square of half stroke, cm2 
short-circuit valve area, cm2 
pivot element in matrix invert routine 
scalar for matrix predictor (biggest 
change) 

scalar for matrix predictor (bi ^gest 
number) 

change of piston 3 position itom initial 
value, cm 


CCCHKP 

CCCKKK 


CMAT 

COSl ^ 

COS2 1 

COS3 I 

COS4 j 

CYCLPM 

CYCLPR 

DEGR 

DELE 

DELPSI 

BELT 

DELTAV 

DELX 

DELY 

DENOM 
DETERM 
DIAM 
DIRl ^ 
DIR2 1 
DIR3 f 
DIR4 J 
DIV 
DRIVO 
E 

EMAT 
ERRBSE 
FDRAG 
FFFFl ^ 
FFFF2 1 
FFFF3 [ 
FFFF4 J 
FORCl ^ 
FORC2 1 
FORC3 f 
FORC4 J 
FRAC 
FREQ 
FREQQ 
FREQRP 
FRRF 


previous change in piston 3 position, cm 
scalar change to time increment to 
predict when piston 3 will again reach 
its initial value 
storage matrix in BROYFl 

cosine of crank angle and 90'* offsets 


maximum cycle pressure, MPa 
current cycle pressure, MPa 
conversion factor, rad/deg 
vector of changes in error vector 
change in crank angle, rad 
change in time, sec 
vector change in guess variables 
vector change in guess variables 
(Broyden algorithm) 
vector change in error variables 
(Broyden algorithm) 
scalar change to cycle pressure, cm3/K 
matrix determinant 
orifice diameter, cm 

piston ring leakage flow direction 


stored mass scalar 

initial condition for drive torque, N-m 
error vector 

Jacobian matrix for 14th-order system 

past error vector 

force due to vehicle air drag, N 

torque scalar for each piston 


force on each piston, N 


external control for matrix convergence 
engine rotational speed, rps 
engine rotational speed, rpm 
engine rotational speed, rpm 
force due to rolling resistance, N 
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G 

gravitational constant, kg/(MPa sec2) 

GR 

gear ratio 

GRAT 

overall gear ratio 

GTRAN 

transmission gear ratio 

HOLD 

variable for row and column 
interchange in matrix invert routine 

I 

integer 

IBRYTH 

switch for using Broyden algorithm 

ICALC 

switch for output options 

ICON 

counter for converged errors 

IDB 

switch for internal DEBUG output 

IHPCNV 

switch for causing a matrix to be 
generated at every time point 

IICC 

counter for partial derivative generation 

IJ 

index 

IK 

index 

INDICA 

switch for sign change on error variable 

IPRTOP 

switch for number of printouts per cycle 

IROT 

switch for indicating a complete crank 
rotation 

ISHSP2 

two time points after short circuiting 
starts 

ISHSTP 

time point for end of short circuiting 

ISHSUP 

one time point after short circuiting 
starts 

ISHTST 

time point for start of short circuiting 

ISPSTP 

time point for end of supply 

ISS 

switch for initial conditions 

ISTEP 

counter 

ISUPST 

time point for start of supply transient 

ITRAN 

counter for time steps 

ITRMAX 
IZ S 

J 
JI 

maximum number of time steps 

JN 

JP 

JQ 

JR J 

^ index 

JWORK 

index for setting time increment for 
crank rotation 

K 

index 

KBROY 
KI -I 

counter for calls to Broyden algorithm 

KJ ^ 

KK J 

index 

KWORK 

counter for pressure-volume integration 

KWRIT 

counter for printout 

L 

index 


LW 

M 

MATRIX 

MATTOT 

MBUG 

MPAS 

MW 

N 

NCYSHS 

NCYSHT 

NCYSTP 

NCYSUP 

NITER 

NK 

NMAT 

NOBUG 

NONENG 

NPSI 

NPTPCC 

NPTPCY 

NSTP 

NTMAX 

NUMBCY 

PAUX 

PAVEMP 

PCCBAR 

PCCMAX 

PCCMIN 

PCMAX 


PCMIN 


PCNCHG 
PCOLDI ^ 
PCOLD2 1 
PCOLD3 f 
PCOLD4 J 


scratch vector 
index 

switch for generating a new Jacobian 
matrix 

counter for number of matrices 
generated during a transient run 
debug switch for iteration failure 
maximum allowable iteration passes 
scratch vector 
system order 

cycle at which to stop short circuiting 
cycle at which to start short circuiting 
cycle at which to stop supply 
cycle at which to start supply 
counter for number of iterations at a 
time point 
index 

counter for number*of Jacobian 
matrices generated at a time point 
switch for a debug printout 
switch for type of piston motion 
reset PSI to zero after 360° rotation 
number of integration points per cycle 
plus 1 

number of integration points per cycle 
switch for storing past converged scale 
factors for iteration guesses 
maximum number of iteration variables 
number of engine cycles to be run 
power lost due to engine auxiliaries, kW 
cycle mean pressure for all four 
cylinders after iteration, MPa 
cycle mean pressure for all four 
cylinders during iteration, MPa 
cycle maximum pressure for all four 
cylinders during iteration, MPa 
cycle minimum pressure for all four 
cylinders during iteration, MPa 
cycle maximum compression-space 
pressure for all four cylinders after 
iteration, MPa 

cycle minimum compression-space 
pressure for all four cylinders after 
iteration, MPa 
iteration convergence rate 

compression-space pressure from 
previous time step for pressure-volume 
integration, MPa 
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PCVCOl 
PCVC02 1 
PCVC03 [ 
PCVC04 J 
PC VCR 1 
PCVCR2 I 
PCVCR3 [ 
PCVCR4 J 
PCY 
PCYl 

PCY2 1 
PCY3 [ 
PCY4 J 
PCI ^ 

PC2 ^ 

PC3 
PC4 ^ 

PCI MAX ^ 
PC2MAX I 
PC3MAX f 
PC4MAX J 
PCIMIN 
PC2MIN I 
PC3MIN f 
PC4MIN J 
PDRAG 
PEMAX 


PEMIN 


PENG 


PEOLDl ^ 
PEOLD2 I 
PEOLD3 ( 
PEOLD4 j 
PEVEOl 
PEVE02 I 
PEVE03 r 
PEVE04 J 



PFRICT 


summation of pressure-volume 
integration in compression space to 
previous time step, kW 

summation of pressure-volume 
integration in compression space to 
current time step, kW 

cycle pressure, MPa 

cycle pressure in each working 
space, MPa 


compression-space pressure in each 
working space, MPa 


maximum compression-space pressure 
in each working space, MPa 


minimum compression-space pressure 
in each working space, MPa 

power loss due to air drag, kW 

cycle maximum expansion-space 
pressure for all four cylinders after 
iteration, MPa 

cycle minimum expansion-space 
pressure for all four cylinders after 
iteration, MPa 

engine power loss due to auxiliaries and 
mechanical friction, kW 

expansion space pressure from previous 
time step for pressure-volume 
integration, M 

summation of pressure-volume 
integration in expansion space to 
previous time step, kW 

summation of pressure-volume 
integration in expansion space to 
current time step, kW 


PIE 

PLOAD 
PNET 
POSDEG 
POSDEO 
POSDT 
POSEOO 
POWERT 
POWERl ^ 
POWER2 I 
POWERS [ 
POWRE4 J 
PPl 

PP2 I 
PP3 [ 

PP4 J 
PRRF 
PRl 1 
PR2 I 
PR3 [ 

PR4 J 
PSHTMN 

PSHTMX 

PSI 

PSIDDT 

PSIDEG 

PSIDT 

PSIM 

PSIMD 
PSIWDT 
PSIl ^ 
PSI2 I 
PSI3 [ 
PSI4 J 

PSOURC 

R 

RAT 

RATIO 

RATOO 


RATT 

expansion-space pressure in each RATTT 

working space, MPa 


power loss due to mechanical engine REF 

friction, kW 


constant (3.14159) 
power loss due to load, kW 
net power, kW 
crank angle, deg 

initial-condition crank angle, deg 
speed, rad/sec 

initial-condition speed, rad/sec 
total engine power, kW 

power from each working space, kW 


average pressure at piston ring, MPa 


power loss due to rolling resistance, kW 

regenerator pressure in each working 
space, MPa 

minimum-pressure manifold 
pressure, MPa 
maximum-pressure manifold 
pressure, MPa 
crank angle, rad 

crank angle acceleration, rad/sec^ 
crank angle, deg 
crank angle velocity, rad/sec 
crank angle adjusted to remain between 
0 and 2x radians, rad 
adjusted crank angle, deg 
wheel velocity, rad/sec 

piston crank position, rad 


working-fluid storage pressure, MPa 
gas constant, (MPa cm3)/(K kg) 
scalar change on step size for iteration 
largest step-size change 
ratio of slot valve opening length to 
piston rod length 
length of slot valve opening, cm 
scalar adjustment to expansion- and 
compression-space temperatures due to 
speed 

desired value of summation of error 
changes 
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RER 

RMAT 

RODL 

RODLS 

RR 

RRC 

RS 

RSHT 

RSHTCT 

RSHTT 

RSUP 

RWHEEL 

SES 

SESP 

SGN 

SlNl 

SIN2 I 
SIN3 [ 
SIN4 J 
SINll S 
SIN22 1 
SIN33 ( 
SIN44 J 

SPEEDR 

SPDMAX 

SPDRPM 

STROKE 

TAUX 

TCOl 

TC02 I 

TC03 f 

TC04 J 

TCI 

TC2 

TC3 

TC4 ^ 

TDRAG 

TEMP 

TEMPI 

TEMP2 

TEMP3 

TENG 


flow resistance between expansion space 
and regenerator, (MPa sec)/kg 
Jacobian matrix for 16th-order system 
crank rod length, cm 
SQRT(RODL**2 + AMPLIT*"2), cm 
maximum error value 
flow resistance betwen regenerator and 
compression space, (MPa sec)/kg 
square of crank length, cm2 
flow resistance between manifolds and 
compression spaces, (MPa sec)/kg 
flow resistance between high- and low- 
pressure manifolds, (cm2 MPa sec)/kg 
short-circuit flow scalar, kg/(sec MPa) 

supply flow resistance, 

(cm2 MPa sec)/kg 

radius of car wheel, cm 
summation of squares of present errors 
summation of squares of past errors 
sign of perturbation step 

square of sines of crank angle and 
offsets 


sine of crank angle and offsets 


engine speed, rad 
maximum engine speed, rpm 
engine speed, rpm 
piston stroke length, cm 
torque loss due to auxiliaries, N-m 

initial value of compression-space 
temperatures for each working 
space, K 


compression-space temperature for each 
working space, K 

torque loss due to air drag, N-m 
Broyden update scalar 
Broyden update scalar 
Broyden update vector 
Broyden update vector 
total torque loss due to engine 
auxiliaries and mechanical friction, 
N-m 


TESTl 
TEST2 I 
TESTS [ 
TEST4 J 
TEOl 
TE02 
TE03 
TE04 

TEl 

TE2 I 
TE3 ( 
TE4 J 

TFRICT 

TIME 
TIMPRV 
TLOAD 
TNET 
TOLPCG 



TOLSS 

TOLl 

TOL2 

TORMAX 

TORMIN 

TORQT 

TORQl 
TORQ2 I 
TORQ3 [ 
TORQ4 J 
TRRF 
TROl ^ 
TR02 I 
TR03 ( 
TR04 J 
TRl 

TR2 I 
TR3 r 
TR4 J 

TSH 

TST31H ^ 
TST32H I 
TST33H ( 
TST34H J 
TST31L S 
TET32L 1 
TST33L r 
TST34L J 


pressure drop across piston rings for 
each piston, MPa 

initial value of expansion-space 
temperatures for each working 
space, K 


expansion-space temperature for each 
working space, K 

torque loss due to mechanical 
friction, N-m 
time, sec 

time at previous time step, sec 
torque loss due to load, N-m 
net torque, N-m 

convergence rate at which a decision is 
made to generate a new Jacobian 
matrix 

error tolerance 

lower limit for good partial derivatives 
upper limit for good partial derivatives 
maximum torque, N-m 
minimum torque, N-m 
total torque, N-m 


torque generated by each piston, N-m 


torque loss due to rolling friction, N-m 

initial-condition regenerator gas 
temperature for all four working 
spaces, K 


regenerator gas temperature for all four 
working spaces, K 


short-circuit manifold temperature, K 


pressure drop from compression spaces 
to high-pressure plenum, MPa 


pressure drop from low-pressure 
plenum to compression spaces, MPa 
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TTl ^ 

TT2 1 
TT3 ( 

TT4 J 
TWC 

TWCIN 

TWH 

TWHIN 
VCCl ^ 
VCC2 1 
VCC3 [ 
VCC4 J 
VCOLDl ^ 
VCOLD2 I 
VCOLD3 [ 
VCOLD4 J 


average temperature at piston ring for 
each piston, K 

adjusted compression-space 
temperature, 

cooler wall temperature, °C 
adjusted expansion-space 
temperature, °C 
heater wall temperature, °C 

initial-condition compression-space 
volume for each working space, cm3 

compression-space volume from 
previous time step for pressure-volume 
integration, cm3 


VMAT 

voc 

VOE 

VR 

VRl ^ 
VR2 [ 
VR3 I 
VR4 J 

VS 

VSAVE 


VSH 

WDSHT 


VCONV 

VCl 

VC2 1 
VC3 f 
VC4 J 
VDELTA 


VDENOM 


VDOT 


VDOTSV 


VDOTT 


VEEl 
VEE2 I 
VEE3 j 
VEE4 J 
VEOLDl ^ 
VEOLD2 I 
VEOLD3 ( 
VEOLD4 J 
VEl ^ 

VE2 I 
VE3 [ 

VE4 J 


vector of converged state variables from 
previous time step 


compression-space volume, cm3 


initial perturbation on state variables 
for matrix generation 
vector of past converged state variables 
used in error vector calculation 
vector of state variable derivatives at 
current time 

vector of state variable derivatives at 
previous time 

vector of average values of state variable 
derivatives from present and past time 
steps 


initial-condition expansion-space 
volume for each working space, cm3 


expansion-space volume from previous 
time step for pressure-volume 
integration, cm3 


expansion-space volume, cm3 


WDl \ 

WD2 / 

WD3 ) 

WD4 J 
WD5 \ 

WD6 J 
WD7 1 
WD8 / 
WDIIN ^ 
WD2IN I 
WD3IN I 
WD4IN J 
WDIOUT ^ 
WD20UT I 
WD30UT f 
WD40UT J 
WLK12 
WLK23 1 
WLK34 r 
WLK41 J 
WSCIDT ^ 
WSC2DT 1 
WSC3DT [ 
WSC4DT j 
WSEIDT ^ 
WSE2DT I 
WSE3DT ( 
WSE4DT J 

WSMNl 


VGUESS state variable guess vector WSMXl 

VKPH car velocity, km/hr 


vector of changes in state variables 
during iteration 

compression-space dead volume, cm3 
expansion-space dead volume, cm3 
regenerator volume, cm3 

regenerator volumes for the four 
working spaces, cm3 

state variable guess vector 
vector of saved, converged state 
variables used for generating a 
Jacobian matrix 

short-circuit low- and high-pressure 
plenum volumes, cm 
short-circuit flow from high-pressure 
plenum to low-pressure plenum, 
kg/sec 

mass flow rates in working space 1 , 
kg/sec 

mass flow rates in working space 4, 
kg/sec 

mass flow rates in working space 3, 
kg/sec 

mass flow rates in working space 2, 
kg/sec 

short-circuit mass flows from low- 
pressure plenum to compression 
spaces, kg/sec 

short-circuit mass flows from 
compression spaces to high-pressure 
plenum, kg/sec 


piston ring leakage flow from working 
space i to j, kg/sec 


stored mass derivative in compression 
volume for each working space, kg/sec 


stored mass derivative in expansion 
volume for each working space, kg/sec 


initial stored mass in low-pressure 
plenum, kg 

initial stored mass in high-pressure 
plenum, kg 
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WSRIDT 
WSR2DT I 
WSR3DT I 
WSR4DT J 

WSSCNl ^ 
WSSCN2 I 
WSSCN3 ( 
WSSCN4 J 
wssn ^ 

WSS22 I 
WSS33 f 
WSS44 j 


stored mass derivative in regenerator 
volume for each working space, kg/sec 


piston rod leakage flow for each 
working space, kg/sec 


stored mass supply flow into each 
working space, kg/sec 


XDST 

XDl 

XD2 I 
XD3 ( 
XD4 J 
XDIP 
XD2P I 
XD3P [ 
XD4P J 
XD30 
XD3PRV 


WSTMN 

WSTMND 

WSTMX 

WSTMXD 

WSTOTL 

WSWANT 

WTCAR 

WTENG 

WTOT 

WTWHEL 


stored mass in low-pressure plenum, kg 
stored mass derivative in low-pressure 
plenum, kg/sec 

stored mass in high-pressure plenum, kg 
stored mass derivative in high-pressure 
plenum, kg/sec 

total stored mass in engine, kg 
initial stored mass in each working 
space, kg 
mass of car, kg 
mass of engine, kg 
total stored mass in engine during 
integration, kg 
mass of wheels, kg 


XXX 

Y 

YOLD 

YSCALE 

YYY 

YYYX 


Z 

ZEROES 


WTl 

WT2 1 stored mass in each working space ZMAT 

WT3 I during integration, kg ZOLD 

WT4 J 


slide valve opening position, cm 

piston position with zero reference at 
midstroke, cm 


piston position with zero reference at 
top of crank, cm 

initial position of piston 3, cm 
position of piston 3 at previous time 
step, cm 

summation of squares of changes in 
errors to maximum error 

present abscissa value for trapezoidal 
integration 

past abscissa value for trapezoidal 
integration 
stored mass scalar 
scalar vector for changes in state 
summation of squares of changes in 
errors to maximum error for debug 
output 

present ordinate value for trapezoidal 
integration 

position bias to change zero reference of 
piston position from top crank to 
center stroke, cm 

Jacobian matrix for 12th-order system 
past ordinate value for trapezoidal 
integration 
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Appendix B 
Analytical Model 

The dynamic equations used to model the four-cylinder controls model were simplified from those 
used to represent the single-working-space model. Both sets of equations are given in reference 2. For 
completeness, the equations of the four-cylinder model are presented here. All equations are in 
FORTRAN. 

Dynamic Equations 

A schematic of the four-cylinder controls model is given in figure 1. The pistons are numbered in the 
order in which they reach top stroke as indicated by the arrows in the piston heads. The dynamic 
equations are 


WSE1DT= - WDl -WLK12 
WSR1DT = WD1-WD2 

WSC1DT = WD2-I-WLK41-I-WSS11-I-WSSCN1+WD1IN-WD1 OUT 
WSE2DT= -WD7-WLK23 
WSR2DT = WD7-WD8 

WSC2DT = WD8 -h WLK12 WSS22 -h WSSCN2 -h WD2IN - WD20UT 
WSE3DT= -WD5-WLK34 
WSR3DT = WD5-WD6 

WSC3DT = WD6 -f- WLK23 -I- WSS33 -I- WSSCN3 -I- WD3IN - WD30UT 
WSE4DT = - WD3 - WLK41 
WSR4DT = WD3-WD4 

WSC4DT = WD4 -h WLK34 -I- WSS44 -I- WSSCN4 -I- WD4IN - WD40UT 


The first-order approximations to the momentum equations are 


WD1=(PE1-PR1)/RER 
WD2 = (PR1 -PC1)/RRC 
WD3 = (PE4-PR4)/RER 
WD4 = (PR4 - PC4)/RRC 
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WD5 = (PE3-PR3)/RER 


WD6 = (PR3 - PC3)/RRC 

WD7 = (PE2 - PR2)/RER 

WD8 = (PR2 - PC2)/RRC 
Pressure 

Pressures are given by the ideal-gas law: 
PEI =(TErR*WSEl)/VEl 
PRl =(TR1*R*WSR1)/VR1 
PCI =(TC1*R*WSC1)/VC1 
PE2 = (TE2*R* WSE2)/VE2 
PR2 = (TR2*R*WSR2)/VR2 
PC2 = (TC2*R*WSC2)/VC2 
PE3 = (TE3*R*WSE3)/VE3 
PR3 = (TR3*R*WSR3)/VR3 
PC3 = (TC3*R*WSC3)/VC3 
PE4 = (TE4*R‘WSE4)/VE4 
PR4 = (TR4*R*WSR4)/VR4 
PC4 = (TC4*R* WSC4)/VC4 


Temperatures 

Temperatures are set as a weak function of speed. The temperature functions in the expansion and 
compression spaces are derived from steady-state data from the single-working-space model results 
(ref. 2). The speed ratio is 


RATTT = (SPDMAX - FREQRP)/1000.0 
and the temperatures in the volumes are given by 
TEl =TWH + RATTT*5.555 - 55.55 
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TE2 =TWH + RATTT*5.555 - 55.55 
TE3 = TWH + RATTT*5.555 - 55.55 
TE4 = TWH + RATTT*5.555 - 55.55 
TCI =TWC-RATTT*3.888 + 55.55 
TC2 =TWC - RATTT*3.888 + 55.55 
TC3 =TWC - RATTT*3.888 + 55.55 
TC4=TWC-RATTT*3.888 + 55.55 
TRl=(TEl+TCl)/2.0 
TR2 = (TE2 + TC2)/2.0 
TR3 = (TE3+TC3)/2.0 
TR4 = (TE4 + TC4)/2.0 
Volumes 

Variable volumes are a function of piston position: 

VE 1 = VOE + (STROKE/2.0 - XD 1 ) * AD 
VCl = VOC + (STROKE/2.0 + XD4)*(AD - AR) 

VE2 = VOE + (STROKE/2.0 - XD2)*AD 
VC2 = VOC + (STROKE/2.0 + XD1)*(AD - AR) 

VE3 = VOE + (STROKE/2.0 - XD3)* AD 
VC3 = VOC + (STROKE/2.0 + XD2)*(AD - AR) 

VE4 = VOE + (STROKE/2.0 - XD4)* AD 
VC4 = VOC + (STROKE/2.0 + XD3)*(AD - AR) 

The regenerator volumes are fixed: 

VR1=VR 
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VR2=VR 


VR3 = VR 

VR4 = VR 
Piston Positions 

Piston positions are calculated from the drive geometry and the crank angle: 

XDl =(STROKE/2.0)*COS(PSI) 

+ SQRT(RODL* *2 - (STROKE/2.0)* *2*(SIN(PSI))* *2) 

- SQRT(RODL* *2 - (STROKE/2.0)* *2) - ZEROES 

XD2 = (STROKE/2.0)*COS(PSI - PIE/2.0) 

+ SQRT(RODL*»2-(STROKE/2.0)**2*(SIN(PSI- PIE/2.0))* *2) 

- SQRT(RODL* *2 - (STROKE/2.0)**2) - ZEROES 

XD3 =(STROKE/2.0)*COS(PSI-PIE) 

+ SQRT(RODL* *2 - (STROKE/2.0)* *2*(SIN(PSI - PIE))* *2) 

- SQRT(RODL* *2 - (STROKE/2.0)* *2) - ZEROES 

XD4 = (STROKE/2.0)*COS(PSI - 3. *PIE/2.) 

+ SQRT(RODL**2-(STROKE/2.0)**2*(SIN(PSI-3.*PIE/2.))**2) 

- SQRT(RODL**2 - (STROKE/2.0)* *2) - ZEROES 

where ZEROES is a small bias to center the piston stroke around zero reference rather than at the top 
of the crank. 

Torque 

Torque is calculated from piston areas, cycle pressures, and drive geometry. The forces on the 
pistons are 

FORCl =PE1 *AD - PC2*(AD - AR) 

FORC2 = PE2* AD - PC3 *(AD - AR) 

FORC3 = PE3*AD - PC4*(AD - AR) 

FORC4 = PE4*AD - PCI *(AD - AR) 

The torques generated by each piston are 
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TORQl =(STROKE/2.0)*FORC1*SIN(PSI)*(1.0 + COS(PSI) 

/SQRT(RODL* *2 - (STROKE/2.0)* »2*(SIN(PSI))* *2)) 

TORQ2 = (STROKE/2.0)*FORC2*SIN(PSI - PIE/2.0)*(1 .0 + COS(PSI - PIE/2.0) 

/SQRT(RODL**2-(STROKE/2.0)**2*(SIN(PSI-PIE/2.0))**2)) 

TORQ3 = (STROKE/2.0)*FORC3 *SIN(PSI - PIE)*(1 .0 + COS(PSI - PIE) 

/SQRT(RODL**2-(STROKE/2.0)**2*(SIN(PSI-PIE))**2)) 

TORQ4 =(STROKE/2.0)*FORC4*SIN(PSI-3.*PIE/2.)*(l. + COS(PSI-3.*PIE/2.) 
/SQRT(RODL* *2 - (STROKE/2.0)* *2‘(SIN(PSI - 3.*PIE/2.))* *2)) 

where total torque is 

TORQT = TORQ 1 + TORQ2 + TORQ3 + TORQ4 
Losses 

Vehicle load effects, engine auxiliary power loss, and mechanical friction are calculated from 

PAUX = 1.57667E - 10*FREQRP**3 -4.3E - 7*FREQRP**2 + 8.9331E 
-4*FREQRP- + .78002 

PFRICT=12.8/20.0*(FREQRP/SPDMAX)*(PCCBAR + 5.0) 
where 

PCCBAR = .462*PCCMAX + .538*PCCMIN 

and PCCMAX and PCCMIN are the average maximum and minimum compression-space pressures 
for each of the four working spaces. 

Rolling resistance and air drag are considered as part of the vehicle load. For rolling resistance it is 
assumed that there is a 4.51-hp loss at 50 mph; thus 

PRRF = (4.5 1 * .7457/(73 .333 *30.48)*GRAT*RWHEEL*PSIDT 
where 

GRAT = GTRAN*(1 .0/GR) 

The force due to air drag is 

FDRAG = .010353*4.448*(RWHEEL*GRAT*PSIDT)**2/(30.48)**2 


The power loss is 



PDRAG = FDRAG*RWHEEL*GRAT*PSIDT/100000.0 


The associated torque losses are 

TFRICT = 33000.0*1.356*PFRICT/(2.0*PIE*.7457*FREQRP) 

TAUX = 33000.0*1 .356*PAUX/(2.0*PIE*.7457*FREQRP) 

TRRF = 4.5 1*550.0* 1 ,356*RWHEEL/(73 .333 *30.48)*GRAT 
TDRAG = FDRAG*RWHEEL*GRAT/100.0 
The vehicle speed is 

VKPH = (3600.0*RWHEEL*GRAT*PSIDT)/100000.0 
The net torque is 

TNET = TORQT - TRRF - TDRAG - TFRICT - TAUX 

From the net torque when the drive dynamics are integrated, the equation of motion is 
AINERT*PSIDDT = TNET 


Supply Flow 

Supply flow into each working space is calculated as a function of differential pressure between the 
compression spaces and the supply. Flow is let in by a timing slot in the piston rod. When the rod is 
near bottom stroke, the slot allows flow into the compression space. A schematic of the system is 
shown in figure 7. The supply check valves function such that 


WSS11 = 


WSS22 = 


WSS33 = 


(PSOURC-PCl)*ADSl/RSUP for PClrsPSOURC 

0.0 for PCl>PSOURC 

(PSOURC-PC2)*ADS2/RSUP for PC2<PSOURC 

0.0 for PC2>PSOURC 

(PSOURC-PC3)*ADS3/RSUP for PC3<PSOURC 


0.0 


for PC3>PSOURC 



Figure 7. - Supply system. 


WSS44 = 


(PSOURC - PC4)*ADS4/RSUP for PC4<PSOURC 

0.0 for PC4> PSOURC 


where 


ADS1 = 


r 1.0- (XD4 + (STROKE/2.0))/(RATOO*STROKE) 
J for XD4 < - STROKE*(.5 - RATOO) 

L 0.0 for XD4>-STROKE*(.5 -RATOO) 


ADS2 = 


1.0 -(XDl +(STROKE/2.0))/(RATOO* STROKE) 
for XDl < - STROKE*(.5 - RATOO) 


0.0 


for XDl > - STROKE*(.5 - RATOO) 














f 1.0- (XD2 + (STROKE/2.0))/(RATOO*STROKE) 


ADS3 = 


ADS4 = 



r 

< 


for XD2 < - STROKE*(.5 - RATOO) 

0.0 for XD2>-STROKE*(.5 -RATOO) 

1 .0 - (XD3 + (STROKE/2.0))/(RATOO*STROKE) 
for XD3 < - STROKE*(.5 - RATOO) 

0.0 for XD3 > - STROKE*(.5 - RATOO) 


Rod leakage flow is also modeled. A schematic of the flow is shown in figure 8. 


WSSCNl = 


ARLEAK*(PSOURC - PC1)/RSUP 


0.0 


for PCl>=PSOURC 


WSSCN2 = 


WSSCN3 = 




ARLEAK*(PSOURC - PC2)/RSUP 

0.0 for PC2> =PSOURC 

ARLEAK*(PSOURC - PC3)/RSUP 


WSSCN4 = 


0.0 for PC3>=PSOURC 

ARLEAK*(PSOURC - PC4)/RSUP 


0.0 


for PC4>=PSOURC 


for PCI sPSOURC 


for PC2<PSOURC 


for PC3SPSOURC 


for PC4<PSOURC 


Piston Rod Leakage 

An orifice equation is used to model piston ring leakage: 

WLK12 = .9*AO*ALEAK*SQRT(2.0*G*PP1/(R*TT1))*SQRT(ABS(TEST1))*DIR1 
WLK23 = .9*AO*ALEAK*SQRT(2.0*G*PP2/(R*TT2))*SQRT(ABS(TEST4))*DIR4 
WLK34 = .9 ‘ AO * ALEAK*SQRT(2.0*G *PP3/(R *TT3)) *SQRT(ABS(TEST3)) *DIR3 
WLK41 = .9*AO‘ALEAK*SQRT(2.0*G*PP4/(R*TT4))*SQRT(ABS(TEST2))*DIR2 
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where 

AO = PIE/4.0*(DIAM* *2 - (DIAM - .0001)* *2) 
and 

DIAM = SQRT(4.0‘ AD/PIE) 

Also 

TEST1=PE1-PC2 
TEST2 = PE4-PC1 
TEST3 = PE3-PC4 
TEST4 = PE2-PC3 
PPl=(PEl+PC2)/2.0 
PP2 = (PE4 + PCl)/2.0 
PP3 = (PE3 + PC4)/2.0 
PP4 = (PE2 + PC3)/2.0 
TTl=(TEl+TC2)/2.0 
TT2 = (TE4+TCl)/2.0 


TT3 = (TE3 + TC4)/2.0 





TT4 = (TE2 + TC3)/2.0 


DIR determines the flow direction: 


DIRI = 


1.0 


for TESTl >0.0 


- 1.0 


for TEST 1< 0.0 


DIR2 = 


1.0 


for TEST2>0.0 


- 1.0 


for TEST2<0.0 


DIR3 = 


1.0 


for TEST3>0.0 


- 1.0 


for TEST3<0.0 


DIR4 = 


1.0 


for TEST4>0.0 


- 1.0 


for TEST4<0.0 


Short Circuiting 

A schematic of a short-circuit system is shown in figure 9. The high 
modeled by constant-temperature volumes. The gas temperature 
compression space. For the check valves 


TST3 1 H = PC 1 - PSHTMX 

TST32H = PC2 - PSHTMX 

TST33H = PC3 - PSHTMX 

TST34H = PC4 - PSHTMX 

TST3 1 L = PSHTMN - PC 1 

TST32L = PSHTMN - PC2 

TST33L = PSHTMN - PC3 

TST34L = PSHTMN - PC4 
The flows are 


- and low-pressure plenums are 
is assumed to be that of the 
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Figure 9. - Short-circuit system. 


WD1IN = 


WD2IN = 


WD3IN = 


WD4IN = 


TST31L/RSHT for TST31L>0.0 

0.0 for TST31L<0.0 

TST32L/RSHT for TST32L>0.0 

0.0 for TST32L<0.0 

TST33L/RSHT for TST33Ls0.0 

0.0 for TST33L<0.0 

TST34L/RSHT for TST34L>0.0 

0.0 for TST34LS0.0 


WD10UT = 


TST31H/RSHT for TST31H>0.0 

0.0 for TST31H<0.0 


WD20UT = 


TST32H/RSHT for TST32HsO.O 

0.0 for TST32H>0.0 
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TST33H/RSHT 


for TST33H>0.0 


WD30UT = 


WD40UT = 


0.0 for TST33H<0.0 

TST34H/RSHT for TST34H>0.0 

0.0 for TST34H<0.0 


The short-circuit valve flow is 


WDSHT = (PSHTMX - PSHTMN)*ASHTT/RSHTCT 
where 

PSHTMX = WSTMX*R*TSH/VSH 
PSHTMN = WSTMN*R*TSH/VSH 


The dynamic equations are 

WSTMXD = WDIOUT + WD20UT -I- WD30UT -I- WD40UT - WDSHT 
WSTMND = - WDIIN - WD2IN - WD3IN - WD4IN -I- WDSHT 



Appendix C 
Flow Charts 


This appendix contains flow charts for the main 
program and all the subroutines in the simulation. Most 
flow charts should be easy to follow, with the possible 
exception of that for subroutine FOURW2. FOURW2 
performs both the incrementing of time and the 
integration of the state equations. The integration scheme 
is implicit. It is a backward-difference integration that 
uses a multivariable Newton-Raphson iteration for 
convergence at a time point. The reason for using this 
type of integration is that it is stable for both large and 
small step sizes. This is important when there a large 
spread in eigenvalues for the system being simulated. 
Although only pressure-flow dynamics are being 
simulated, the addition of temperature dynamics may be 
desirable (ref. 2). This would add slow dynamics to the 
simulation. In that event, the implicit integration scheme 
will be available to handle the widespread dynamics while 
insuring stability. 

To help in understanding how subroutine FOURW2 
works, the following description is provided. Statement 
numbers corresponding to the FORTRAN listing are 
given in the flow chart. 


Integration Scheme 

The integration scheme is a backward-difference 
method that uses a multivariable Newton-Raphson 
iteration scheme for convergence. State variables are 
updated by using the old state vector YSoy, the current 
error vector E, and the inverse of a Jacobian matrix 
EMAT (ref. 5). 

VSnew= -EMAT-lxE-l-VSoid (Cl) 

EMAT is a Jacobian matrix of partial derivatives (i.e., 
changes in error variables with respect to changes in state 
variables): 

EMAT(I, J) = (E(J) - ERRBSE(J))/DELTAV(I) (C2) 

Updating takes place when the errors are converged 
within tolerance. With this technique, both steady-state 
and transient solutions can be obtained by changing the 
error variables. In steady state, all states are at rest; thus 

VDOT = 0.0 = E (C3) 

For a transient case 


VDOTxDELT-(VSnew-VSoid) = E (C4) 

Equations (C3) and (C4) are converged when all the 
elements of E are within a specified tolerance. 

To use this method, a Jacobian matrix must be 
calculated (usually by finite differences) and then 
inverted. This is usually very time consuming. Therefore 
the logic in FOURW2 allows for calculation of a new 
matrix only under adverse conditions such as when the 
rate of convergence of the simulation is getting too slow 
(PCNCHG less than TOLPCG), or when the number of 
allowable passes (MPAS) has been exceeded. 

Perturbation Calculation 

There are several features in FOURW2 that help the 
implicit integration scheme converge. Of primary 
importance is the generation of a “good” Jacobian 
matrix. All the partial derivatives must be representative 
of the linear behavior of the system at a given operating 
point. Since finite differences are used, the sizes of the 
perturbations of the states are important. If they are too 
large, errors will be introduced by the system 
nonlinearities. If they are too small, the partials will be in 
error because of numerical problems (without double- 
precision arithmetic). Thus a tuning mechanism has been 
included in FOURW2 to optimize the sizes of the 
perturbations. First, the sum of squares of all the changes 
in the errors is calculated for each perturbation. Once this 
is done, the “goodness” of the partial is checked by 
calculating for each state variable 

XXX = r,V E [(E(I)-ERRBSE(I)]2 (C5) 

A > i 

and then checking if 

TOL 1 < XXX <TOL2 (C6) 

If all XXX’s fall within the tolerance band, the matrix is 
considered “good”. For this simulation, TOLl =0.001 
and TOL2 = 0.01. Since all the errors are scaled, this 
tolerance band lies between 0.1 and 1 percent. For a more 
linear system the band could be larger; for a more 
nonlinear system, smaller. 

Scaling of Perturbations 

In general, for the initial perturbations at a point, the 
XXX’s will not fall within the tolerance band described 
above. Thus, FOURW2 scales the perturbations to try to 
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force the XXX’s within the band. This is done by 
calculating 


YYY = REF/XXX 


(C7) 


for each state variable. REF is defined as being in the 
center of the tolerance band: 


REF = (TOLI +TOL2)/2.0 


(C8) 


Once a set of YYY’s has been calculated such that the 
XXX’s fall within the band, the set of YYY’s is stored. 
After this has been done for all N states, the scaling 
vector YYY is generated. When a new matrix is needed, 
the scaling vector YYY is applied to the current states to 
determine first guesses for the perturbations needed to 
obtain new partial derivatives. If for any state variable 
the new XXX falls outside the tolerance band, YYY is 
updated and the new result stored. This method generally 
reduces the number of passes required for subsequent 
matrix generation. 

Error Messages 

In generating a partial derivative a situation may arise 
where XXX never gets within tolerance. When this 
happens, the program prints out an error message 
“CHECK INPUT -BAD PARTIAL DERIVATIVE,” 
prints out a debug output to help the user diagnose the 
problem, and then stops the simulation. This is the only 
time when the simulation is stopped except for the 
normal exit when ITRAN has been incremented to its 
maximum value (ITRMAX). In general, this problem will 
occur when coding is added to the simulation that is not 
consistent throughout the simulation. One example 
would be calculating a piston ring leakage flow from one 
working space but neglecting to add it to the adjacent 
working space. 

Broyden Update Algorithm 

As stated earlier, it is time consuming to calculate a 
Jacobian matrix and even more time consuming to invert 
the matrix for a large system. Usually double precision is 
required (as in DMINV). Lack of speed in the inverse 
algorithm can be prohibitive, especially in a controls 
model, where long transients are usually needed to 
evaluate controls schemes. One possible means of 
avoiding this problem is to generate the inverse and then 
to update it continually with information gained as the 
inverse is used. This method is the Broyden update 
algorithm (ref. 6) and is found in subroutine BROYFl. 

However, it should be noted that using the method is 
also time consuming and, for low-order systems, may not 


help much in reducing the computation time for the 
simulation even when using large time steps. Table XIX 
contains the matrix update strategy for the simulation. If 
the Broyden algorithm is not used, a new Jacobian is 
generated only when the convergence rate (PCNCHG) 
falls below 0.5. At values of PCNCHG above 0.5 the 
original inverse is used. When using the Broyden 
algorithm a new Jacobian matrix is calculated only when 
TOLPCG falls below 0.0 or when the algorithm has been 
called more than twice the order of the system plus one 
(2*N+1) times. Since it is time consuming to use the 
algorithm, the inverse is updated only when PCNCHG 
falls below 0.7. Above PCNCHG equal to 0.7 the inverse 
is not updated. 

Even though the algorithm does not shorten computer 
time in the current simulation, it has been included in the 
program. If the energy equation is added to the volume 
dynamics or if a control system is added such that the 
number of states in the simulation increases significantly, 
the algorithm will prove much more useful. 

Main program MAINSE 


Input engine 
geometry and healer 
and cooler wall 
temperatures 


Input flow 
resistances, matrix 
convergence data, and 
switch settings 


Input run conditions, 
load characteristics, 
cycle data, and 
transient desired 


Calculate number of 
lime steps, supply transient 
start, short-circuit start, 
and inertias 
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Subroutine BROYFl 



Subroutine DMINV 









































Check to see if perturbation 
size caused significant change 
in error vector; calculate 























































Subroutine FWS3V2 



Calculate flow 
and crank angle 
derivatives 


: 

Set VDOl 

o 

o* 

II 



Calculatf 
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initial- 

errors 
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TOURW2 















Subroutine GUESF2 


Subroutine ICSTUP 




FOURW2 




















Subroutine PISTN3 



Subroutine PLEAK2 



43 







Subroutine SHORT2 


















Subroutine SUPPL2 



Subroutine TRAP 




TOURW2 














Appendix D 
Steady-State Results 

Steady-state results using the simulation are presented 
in figures 10 to 13. These results are presented here to 
help the user understand the mode of program operation. 



Figure 10, - Net power as a function of 
speed for single- and four-working- 
space models and engine data. 



Figure 11. - Net power as a function of 
speed for various piston ring leakage 
areas. Engine mean pressure, 
PAVEMP, 15 MPa. 
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Figure 12. - Net power as a function of 
engine speed for various short-circuit 
valve areas. Engine mean pressure, 
PAVEMP, 15 MPa. 
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Figure 13. - Net torque as a function of 
engine speed for various short-circuit 
valve areas. Engine mean pressure, 
PAVEMP, 15 MPa. 



The term “steady state” for a Stirling engine actually 
refers to a transient run in which the engine is simulated 
over many cycles at constant input conditions (constant 
heater and cooler wall temperatures and constant speed; 
i.e., no supply fluid). Figures 10 to 13 were also presented 
in reference 2. All steady-state cases were run by forcing 
piston position as a function of time (NONENG = 0). 

Figure 10 shows a comparison of results from the 
simplified four-cylinder Stirling engine controls model 
with results from the seven-volume, single-working-space 
model (which includes the energy equation). Also shown 
are experimental data run at Lewis. The simulation data 
were obtained by setting constant engine mean pressure 
(CYCLPR = 5., 10., or 15. MPa) and constant engine 
speed (SPDRPM = 1000., 2000., 3000., or 4000. rpm). 
After steady state was reached, cycle data were recorded. 
Figure 10 shows that the simplified model agrees well 
with the experimental data over the ranges of engine 
speed and mean pressure. 


Steady-state power versus engine speed at different 
piston ring leakage areas is given in figure 11. To obtain 
these data, mean pressure was held constant at 15 MPa 
(CYCLPR = 15.0), while speed and leakage area were 
varied (SPDRPM = 1000., 2000., 3000., and 4000. and 
ALEAK = 0.0, 0.002, 0.004, 0.006, 0.008, and 0.010, 
respectively). As expected, performance degrades as 
leakage area increases. However, power loss due to 
leakage seems to be independent of speed. 

Net power versus speed for different short-circuit valve 
areas is shown in figure 12. These curves were generated 
by setting engine mean pressure constant 
(CYCLPR= 15.0) and varying engine speed (SPDRPM) 
and short-circuit valve area (ASHTT). Note the large 
drop in engine power with short circuiting. The 
corresponding torque curve is given in figure 13. Torque 
reversal occurs for leakage areas above 0.1. 
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Appendix E 
Transient Results 


For controls analysis, transient data over a number of 
engine cycles are of interest. To obtain these data, the 
simulation calculates each cycle individually. Thus the 
simulation can also be used to obtain information 
regarding engine parameters during a single cycle. Data 
on an individual cycle basis were presented in reference 4. 
Data are presented herein to help the user understand 
how to obtain transient information. 

A working-fluid supply transient for six cycles is shown 
in figure 14. Because for this case speed must change, the 
simulation must be run with the drive dynamics 
(NONENG = l). The initial engine pressure is 5 MPa 
(CYCLPR = 5.0); the source pressure for the working 
fluid is 13.8 MPa (PSOURC = 13.8); the supply is started 
at the second engine cycle (NCYSUP = 2); there is no rod 



Figure 14. - Working-fluid supply transient with no piston 
rod leakage. Working-fluid storage pressure, PSOURC, 
13. 8 MPa at 0.06 sec. 


leakage (ARLEAK = 0.0). To get the desired printout, set 
IPRTOP equal to 1; this gives a printout at every point 
during the cycle (NPTPCY = 100). 

The same transient with piston rod leakage is shown in 
figure 15. In this case, all the input data are the same 
except ARLEAK =0.5. Note that the transient results for 
the two cases are markedly different. 




E S 
o o. 


o O E 

Qi I 

15 P ^ 




Figure 15. - Working-fluid supply transient with uniform 
piston rod leakage. Piston rod leakage area, ARLEAK, 
0.5 cm^; working-fluid storage pressure, PSOURC, 

13. 8 MPa at 0.06 sec. 
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TABLE I. - SIMULATION INPUT: ENGINE GEOMETRY 


Name 

Setting 

Function 

AD 

23.7613 

Piston area 

AR 

1.1290 

Piston rod area 

RODL 

10.0 

Crank rod length 

VR 

115.339 

Regenerator volume 

VOE 

63.87 

Expansion-space dead volume 

VOC 

61.70 

Compression-space dead volume 

STROKE 

4.0 

Piston stroke 

ALPHA 

90.0 

Crank angle lag 


TABLE II. - SIMULATION INPUT: HEATER 
AND COOLER WALL TENPERATURES 


Name 

Setting 

Function 

TWHIN 

TWCIN 

705.0 

86.0 

Heater wall temperature 
Cooler wall temperature 


TABLE III. - SIMULATION INPUT: 
FLOW RESISTANCES 


Name 

Setting 

Function 

RER 

1.9 

Resistance between expansion 



space and regenerator 

RRC 

1.9 

Resistance between regenerator 

j 


and compression space 


TABLE IV. - SIMULATION INPUT: CONSTANTS 


Name 

Setting 

Function 

DEGR 

PIE 

R 

G 

57.296 

3.1416 

4125.6 

10017.0 

Degrees to radians 

Gas constant 
Gravitational constant 







TABLE V. - SIMULATION INPUT: MATRIX CONVERGENCE 


Name 

Setting 

Function 

VDELTA 

0.01 

Initial perturbation of guesses, 
1 percent 

FRAC 

1.0 

External control of iteration 
step magnitude 

T0L1 

0.001 

Bottom limit on error tolerance 
for matrix linearity 

T0L2 

0.01 

Top limit on error tolerance for 
matrix linearity 

TOLSS 

0.0001 

Solution tolerance 

N 

12 

System order when NONENG = 0 


14 

System order when NONENG = 1 

NTMAX 

16 

Largest system order (when short 
circuiting) 

MPAS 

20 

Maximum allowable convergence 
passes 

TOLPCG 

0.5 

Switch for calculating new matrix 

TABLE VI. - 

SIMULATION INPUT: SWITCHES 

Name 

Setting 

Function 

ISS 

0 

Set up initial conditions 


1 

Transient (set internally) 

ICALC 

0 

Detailed printout 


1 

Short printout 

MATRIX 

1 

Generate a new Jacobian matrix 

NONENG 

0 

Force crank angle as a function 
of time 


1 

Runs with drive dynamics in 

IPRTOP 

NN 

Print out data every NN points 

IBRYTH 

0 

Do not use Broyden algorithm 


1 

Use Broyden algorithm 

IHPCNV 

0 

Use logic to generate a new 
matrix 


1 

Generate a new matrix at every 
point 

NOBUG 

0 

No debug output^ 


1 

Desire debug output 


^If iteration convergence fails or matrix 
generation problems occur, a debug output 
is given. 
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TABLE VII- - SIMULATION INPUT: RUN CONDITIONS 


Name 

Setting 

Function 

ALEAK 

0.0 

Piston ring leakage area scalar 

ARLEAK 

0.0 

Piston rod leakage area 

PSOURC 

10.0 

Hydrogen bottle pressure 

CYCLPM 

15.0 

Maximum engine pressure 

SPDMAX 

4000.0 

Maximum engine speed 

CYCLPR 

5.0 

Initial engine pressure 

POSDEO 

270.0 

Initial crank angle 

WSTOTL 

0.007786 

Total stored mass at 15 MPa 

VSH 

32.78 

Short-circuit volume 

RSHT 

38.0 

Flow resistance between plenums 
and compression spaces 

ASHTT 

0.2 

Short-circuit valve area 


TABLE IX. - SIMULATION INPUT: CYCLE DATA 


Name 

Setting 

Function 

NUMBCY 

NPTPCY 

100 

100 

Desired number of cycles 
Integration points per cycle 


TABLE VIII. - SIMULATION INPUT: 
LOAD CHARACTERISTICS 


[Vehicle information values given in this 
report are representative and do not 
correspond to any actual vehicle-] 


Name 

Setting 

Function 

GTRAN 

GR 

WTENG 

WTWHEL 

WTCAR 

RWHEEL 

1.0/2.53 

1.5 

0.0 

20.48 
1420.0 

30.48 

Transmission gear ratio 

Gear ratio 

Engine mass 

Mass of a wheel 

Car mass 

Wheel radius 


TABLE X- - SIMULATION INPUT: TRANSIENT DESIRED 


Name 

Setting 

Function 

NCYSHT 

5 

Start of short circuit at fifth 



cycle 

NCYSUP 

^15000 

Start of supply at 15 000th cycle 

NCYSTP 

^40000 

End of supply at 40 000th cycle 

NCYSHS 

500 

End of short circuiting at 500th 



cycle 


^Note that if no supply (or short-circuiting) 
is desired, set the corresponding start and 
end parameters very large (greater than the 
number of cycles desired). 






TABLE XI. - COMMON BLOCK LOCATIONS 


Common 

block 

Subroutine 


MAINSE 

BROYFl 

DMINV 

FOURVC 

FWS3V2 

GUESF2 

ICSTUP 

LOSSES 

OUTPUT 

PISTN3 

PLEAK2 

SH0RT2 

SUPLRK 

SUPPL2 

TRAP 

ANGLS 




B 

B 


B 

m 

X 

X 



X 

X 


CARLOD 

X 



WM 

mm 


WM 









CONST 

X 



■■ 

■■ 


■■ 




X 

X 




CYDATA 

X 



WM 

■■ 


WM 


X 







DERIVT 





X 


mM 


X 







ITCONV 


X 


X 

X 




X 

X 


X 




MATRX 


X 


X 




■■ 








MATXIN 

X 

X 


X 

X 


X 









OUTPT 

X 



X 

X 


X 

X 




X 

X 

X 


PCBARR 




X 

X 


X 









PLEKS 





X 


X 




X 





PST 

X 



X 

X 


mm 


X 

X 

X 



X 


RUNCON 

X 



X 

X 


WM 

X 

X 


X 


X 

X 


SAVIT 




X 

X 

X 



X 







SHRTCT 

X 



X 

X 


X 


X 







STEP IT 




mm 

X 


X 

X 








STRIC 




WM 

X 


mM 


X 



^B 

X 

X 


SUPLY 




WM 



WM 


X 




X 

X 


SWITCH 

X 

X 


■■ 

mm 


^B 

X 

X 







TEMPS 





^B 


WM 


X 






V 

TORSUM 




X 



mm 

X 

X 


Hj^B 

^B 




TRANDS 

X 



mm 



■b 






X 

X 


VOLS 




B 

■ll 

IH 

B 


X 


IBI 

■ 







TABLE XII. - 

SAMPLE OUTPUT FOR SHORT 

PRINTOUT 

OPTION 

(I CALC = 

1) 


RUK 

CONDITIONS 

FOR THIS TR&MSIEHT 








TUH = 
ALEAK 

97S.0 
= 0.0000 

TUC = 
ARLEAK 

359.0 
= 0.0000 

CYCLPR = 
PSOURC = 

5.000 

10.00 

NONENG 

ISUPST 

= 

1 

251 

NUMBCY - 
ISHTST = 

1000 

2500001 

TIME 

POSDEG 

POWERT 

TORQT 

PFRICT 

TFRICT 

PAUX 

TAUX 

PRRF 

TRRF 

PDRAG 

TDRAG 

PLOAD 

TLOAD 

PENG 

TENG 

PNET 

TNET 

ITRAN 

KUORK 

PAVEMP 

SPDRPM VKPH 


























TABLE XIII. - SAMPLE OUTPUT FOR DETAILED PRINTOUT OPTION (ICALC = 0) 


STIRLING ENGINE FOUR CYLINDER, THREE VOLUMES PER CYLINDER, CONTROLS MODEL 

1 

RUN CONDITIONS FOR THIS TRANSIENT 


TWH = 

978.0 

TWC = 359.0 

CYCLPR = 

5.000 

NONEHG : 

= 

1 NUMBCY = 

1 

ALEAK = 

0.0000 

ARLEAK = ( 

I.OOOO 

PSOURC = 

10.00 

ISUPST : 

= 501 ISHTST = 5000001 

TIME 

XDl 

XD2 

XD3 

XD4 

POSDEG 

SPDRPM 

PAVEMP 

WTOT DELT 

ITRAN 


PEI 

PRl 

PCI 

TEl 

TRl 

TCI 

VEl 

VRl VCl 

UTl 


PE2 

PR2 

PC2 

TE2 

TR2 

TC2 

VE2 

VR2 VC2 

14T2 


PE5 

PR3 

PC3 

TE3 

TR3 

TC3 

VE3 

VR3 VC3 

WT3 


PE^ 

PR4 

PC4 

TE4 

TR4 

TC4 

VE4 

VR4 VC4 

UT4 


MDl 

UD2 

W55H 

WSSCNl 

14D1IN 

WDIOUT 

WLK12 

WSEIDT WSRIDT 

WSCIDT 


WD7 

UD8 

WSS22 

US5CN2 

WD2IN 

WD20UT 

WLK23 

WSE2DT WSR2DT 

WSC2DT 


WD5 

UDS 

WS533 

W55CN3 

WD3IN 

WD30UT 

WLK34 

WSE3DT W5R3DT 

W5C3DT 


WD3 

UD4 

W5S44 

WSSCN4 

WD4IN 

WD40UT 

WLK41 

WSE4DT WSR4DT 

WSC4DT 


VS(1) 

V5(2) 

V5(3) 

V5(13) 

Ed) 

E(2) 

E(3) 

E(13) TORQl 

POWERl 


VS(^) 

VS(5) 

V5(6) 

VS(14) 

E(4) 

E(5) 

E(6) 

E(14) T0RQ2 

P0WER2 


VS(7) 

VS(8) 

VS(9) 

V5(15) 

E(7) 

E(8) 

E(9) 

E(15) T0RQ3 

P0WER3 


VS(IO) 

VS(ll) 

V5(12) 

VS(16) 

EdO) 

Edl) 

E(12) 

E(16) T0RQ4 

P0WER4 


MSTMXD 

W5TMND 

WD5HT 

PSHTMX 

P5HTMN 

T5H 

WSMXl 

WSMNl 



POWERT 

PFRICT 

PAUX 

PRRF 

PDRAG 

PLOAD 

PENG 

PNET 


POSDEG 

TORQT 

TFRICT 

TAUX 

TRRF 

TDRAG 

TLOAD 

TENG 

TNET KWORK 

SPDRPM VKPH 

0.3000E- 

-01-0.2031 

-2.000 

-0.2008 

2.000 

270.0 

2000. 

5,003 

0.2609E-02 0.3000E-03 

101 


A. 031 

4.015 

4.004 

933.6 

670.2 

406.8 

116.2 

115.3 152.2 

0.6523E-03 


4.552 

4.539 

4.511 

933.6 

670,2 

406 .8 

158.9 

115.3 102.4 

0.6523E-03 


5.985 

6.007 

6.017 

933.6 

670.2 

406.8 

116.2 

115.3 61.70 

0.6523E-03 


5.445 

5.453 

5.481 

933.6 

670.2 

406.8 

63.87 

115.3 102.4 

0.6523E-03 


0.8571E-02 

0.59S7E-02 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-0.8571E-02 0.2584E-02 

0.5987E-02 


0.7222E-02 

0.1462E-01 

0.0000 

0 . 0000 

0.0000 

0 .0000 

0.0000 

-0.7222E-02-0.7403E-02 

0.1462E-01 


-0.1133E-01 

-0.5490E-02 

0.0000 

0 .0000 

0.0000 

0 .0000 

0.0000 

0.1133E-01-0.5836E-02 

-0.5490E-02 


-0.4178E-02 

-0.1491E-01 

0.0000 

0 . 0000 

0.0000 

0.0000 

0.0000 

0.4178E-02 0.1073E-01 

-0.1491E-01 


0.1216E-03 

0.1675E-03 

0.3632E-03 

4.712 

0.9936E-05 

-0.1462E-05 

-0.1090E-05 

0.9294E-07 12.60 

2.671 


0.1878E-03 

0.1893E-03 

0.2752E-03 

209.4 

0.1049E-05 

-0.6665E-05 

0.2206E-05 

0.4292E-07-0.2490E-01 

2.664 


0.1805E-03 

0.2506E-03 

0.2212E-03 

0.1172E-03 

-0.8438E-04 

0.2859E-04 

-0.1996E-04 

0.0000 36.32 

2.662 


0.9029E-04 

0.2275E-03 

0.3345E-03 

0.7813E-04 

0.1733E-05 

-0.1554E-04 

0.6429E-05 

0.0000 -0.5058E-01 

2.671 


0.0000 

0.0000 

0.0000 

6.000 

4.000 

0.0000 

0.1172E-03 

0.7813E-04 



10.67 

3.221 

2.108 

2.531 

2.360 

4.891 

5.329 

0.3469 


270,0 

48.84 

15.38 

10.07 

12.09 

11.27 

23.35 

25.45 

0.4489E-01 101 2000.06 60.56 



LA 

o^ 


TABLE XIV. - SAMPLE DEBUG PRINTOUT (NOBUG = 0) 


- 0.20196623 


- 1.9998884 

- 0.20191854 

2.0001106 






12.062996 


12.062994 

12.062994 

0.00000000 


0.00000000 

116.19153 

152.23161 


13.569489 


13.569489 

13.569489 

0.00000000 


0.00000000 

158.91251 

102.39359 


17.951859 


17.951859 

17.951859 

0.00000000 


0.00000000 

116.19040 

61.702515 


16.429947 


16.429947 

16.429947 

0.00000000 


0.00000000 

63.867355 

102.39467 


922.44482 


668.49976 

414.55493 

922.44482 


668.49976 

414.55493 



922.44482 


668.49976 

414.55493 

922.44482 


668.49976 

414.55493 



0.00000000 


0.00000000 

0.00000000 

0.00000000 


0.00000000 

0.00000000 



0.00000000 


0.00000000 

0.00000000 

0.00000000 


0.00000000 

0.00000000 

0.00000000 

0.00000000 

0.00000000 


0.00000000 

0.00000000 

0.00000000 


0.00000000 

0.00000000 

0.00000000 

0.00000000 

16.000000 


14.000000 

0.00000000 

0.00000000 


0.00000000 




DEBUG OUTPUT 

FRON F 0 URW 2 NUHBER 1 







TIHC ■ 

0.00000000 








DEBUG PRINTOUT 1 

FROH FWS 3 V 2 NUHBER 1 







- 0.20196623 


- 1.9998884 

- 0.20191854 

2.0001106 






12.062987 


12.062986 

12.062989 

0 .S 0193393 E- 

-06 

- 0 . 15058013 E -05 

116.19153 

152.23161 


13.569483 


13.569486 

13.569482 

0.00000000 


0.00000000 

158.91251 

102 . 3935 V 


17.951843 


17.951843 

17.951843 

0.00000000 


0.00000000 

116.19040 

61.702515 


16.429932 


16.429932 

16.429932 

- 0 . 15058013 E- 

-05 

0 . 20077350 E -05 

63.867355 

102.39467 


922.44482 


668.49976 

414.55493 

922.44482 


668.49976 

414.55493 



922.44482 


668.49976 

414.55493 

922.44482 


668.49976 

414.55493 



41 B.B 7988 


7.9077196 

80.492523 

150.37494 


21.255264 

2.6879129 



0.00000000 


0.00000000 

0.00000000 

0.00000000 


0.00000000 

0.00000000 

0.00000000 

0.00000000 

0.00000000 


0.00000000 

0.00000000 

0.00000000 


0.00000000 

0.00000000 

0.00000000 

0.00000000 

16.000000 


14.000000 

0.00000000 

0.00000000 


0.00000000 




DEBUG PRINTOUT FROM FUS3V2 NUHBER 2 








1 

0 . 36830036 E- 03 ' 

- 0 . 81770270 E -06 

0.00000000 


0.00000000 

0.00000000 

0.00000000 

0.00000000 


2 

0 .S 04478 B 1 E -03 

0 . 23878920 E -05 

0.00000000 


0.00000000 

0.00000000 

0.00000000 

0.00000000 


3 

0 . 10737202 E- 02 - 

• 0 . 84144904 E -06 

0.00000000 


0.00000000 

0.00000000 

0.00000000 

0.00000000 


4 

0 . 56662294 E -03 

0 . 1594 S 006 E-OS 

0.00000000 


0.00000000 

0.00000000 

0.00000000 

0.00000000 


5 

0 . 567481 15 E- 03 - 0 . 3714 B 757 E -05 

0.00000000 


0.00000000 

0.00000000 

0.00000000 

0.00000000 


6 

0 . 81239 S 24 E -03 

0 . 14 B 2826 SE -05 

0.00000000 


0.00000000 

0.00000000 

o.oooooooo 

0.00000000 


7 

0 . 54 B 09055 E -03 

0.00000000 

0.00000000 


0.0 '000000 

0.00000000 

0.00000000 

0.00000000 


8 

0 . 75075356 E -03 

0.00000000 

0.00000000 


0.00000000 

0.00000000 

0.00000000 

0.00000000 


9 

0 . 64765452 E -03 

0.00000000 

0.00000000 


0.00000000 

0.00000000 

0.00000000 

0.00000000 


10 

0 . 27573225 E -03 

0.00000000 

0.00000000 


0.00000000 

0.00000000 

0.00000000 

0.00000000 


11 

0 . 687106 S 1 E -03 

0.00000000 

0.00000000 


0.00000000 

0.00000000 

0.00000000 

0.00000000 


12 

0 . 9836 S 918 E -03 

0.00000000 

0.00000000 


0.00000000 

0.00000000 

0.00000000 

0.00000000 


13 

4.7123709 

0.00000000 

0.00000000 


0.00000000 

0.00000000 

0.00000000 

0.00000000 


14 

418.87988 

0.00000000 

0.00000000 


0.00000000 

0.00000000 

0.00000000 

0.00000000 


CONVEROCNCe IN ERR VECTOR 



TABLE 


100 

200 

300 

^00 

500 

600 

700 

800 

900 

1000 

1100 

1200 

1300 

1^00 

1500 

1600 

1700 

1800 

1900 

2000 

2100 

2200 

2300 

2^00 

2500 

2600 

2700 

2800 

2900 

3000 

3100 

3200 

3300 

3^00 

3500 

3600 

3700 

3800 

3900 

^000 

<*100 

<^200 

<♦300 

<*<*00 

<♦500 

<♦600 

<*700 

<*800 

<*900 

5000 

5100 

5200 

5300 

5<*00 

5500 

5600 

5700 

5800 

5900 

6000 

6100 

6200 

6300 

6<t00 

6500 

6600 

6700 


XV. - INPUT ROUTINE (MAINSE) FOR SUPPLY TRANSIENT (100 POINTS/CYCLE) 


INPUT DATA FOR THE FOUR WORKING SPACE THREE VOLUME 
STIRLING ENGINE MODEL. 

COMMON/PST/AD,AR,RODL#VOE»VOC,VR,RER,RRC, STROKE, ALPHA, PSIMD 

COMMON/CARLOD/GTRAN,GR,AIE,AIWHEL,AIVEH,RWHEEL,AINERT 

COMMON/COHST/R,G.PIE,DEGR 

COMMON/MATXIN/VDELTA, FRAC,TOL1,TOL2,TOLSS,N,NTMAX,MPAS, 

1 TOLPCG,REF 

COMMON/SWITCH/ISS.ICALC,MATRIX,NONENG,IPRTOP,IBRYTH. 

1 IHPCNV, NOBUG 

COMMON /SHRTCT/ WDIIN . WD2IN , WD3IN ,WD<*IN ,WD10UT ,WD2OUT,WD30UT,WD<*0U 
1T,WDSHT,WSTMXD,WSTMND,PSHTMX,P5HTMN,WSMX1 ,WSMN1,VSH,RSHT,RSHTT 
COMMON/RUNCON/ALEAK.ARLEAK,PSOURC,CYCLPR,SPDRPM,POSDEO, 

1 WST0TL,TWHIN,TWCIN,CYCLPM,5PDMAX 

COMMON/CYDATA/NPTPCY, ITRMAX , NUMBCY 

COMMOH/TRANDS/ISUPST.ISHTST,TIME,KWORK,ISPSTP,ISHSTP 
COMMON/OUTPT/POSDEG.WTOT, ITRAN,WT1,WT2.WT3.WT<*,WD1 ,WD2,WD3,WD<i, 
1WD5,WD6,WD7 ,UD8,T0RQ1 ,TORQ2,TORQ3,TORQ<*,POWER1,POWER2, POWERS, 
2P0WER<*,P0WERT,T0RQT,TL0AD,TENG,PL0AD,PENG,PNET,PAVEMP,CCCHK 


ENGINE GEOMETRY 
AD=23.7613 
AR=1.1290 
RODL=10.0 
VR=115.339 
VOE=63,87 
VOC=61.70 
STROKE = <*.00 
ALPHA=90.0 

HEATER AND COOLER WALL TEMPERATURES 
TWHIN=705.0 
TWCIN=86.0 

FLOW RESISTANCES BETWEEN VOLUMES 
RER=1.9 
RRC=1.9 

CONSTANTS 
DEGR=57.296 
PIE = 3.1<*16 
R=<*125.6 
G=10017.0 

MATRIX CONVERGENCE INPUT 
VDELTA=.01 
FRAC=1.0 
TQL1=.001 
TOL2=.01 
TOLSS=.0001 
N = l<* 

NTMAX=16 

MPAS=20 
TOLPCG=.5 

SWITCHES 
ISS = 0 
ICALC=1 
MATRIX=1 
NONENG=l 
IPRTOP=100 
IBRYTH=0 
IHPCNV=0 
NOBUG=l 
C 

RUN CONDITIONS 



6fiOQ 


ALEAK=0.0 

6900 


ARLEAK=0.0 

7000 


PSOURC=10.0 

7100 


CYCLPt1=15.0 

7200 


SPDMAX=9000.0 

7300 


CYCLPR=5.0 

7^00 


SPDRPM=2000. 0 

7500 


POSDEO=270.0 

7600 


WSTOTL=. 007786 

7700 


VSH=32.78 

7800 


R5HT=38.0 

7900 


ASHTT=.20 

8000 

C 


8100 

C 

LOAD CHARACTERISTICS 

8200 


GTRAN=1.0/2.53 

8300 


GR=1.5 

8^00 


WTENG=0.0 

8500 


WTWHEL=29.98 

8600 


WTCAR=1920.0 

8700 


RWHEEL=30.98 

8800 

C 


8900 

C 

CYCLE DATA 

9000 


NPTPCY=100 

9100 


NUMBCY=1000 

9200 

C 


9300 

C 

TRANSIENT DESIRED 

9<«00 


NCYSHT=50000 

9500 


NCYSUP=5 

9600 


NCYSTP=50000 

9700 


NCYSHS=50000 

9800 

c 


9900 

c 

CALCULATED INPUT 

10000 


PSHTMX = CYCLPR+1. 0 

10100 


PSHTMN=CYCLPR-1. 0 

10200 


IF (NONENG .EQ. 0) N=12 

10300 


REF=(TOL1+TOL2)/2.0 

10900 


RUHELn=RUHEEL/100.0 

10500 


ITRMAX=NPTPCYXNUMBCY+1 

10600 


ISUPST=NCYSUP«NPTPCY+1 

10700 


ISHTST=NCYSHT«NPTPCY+1 

10800 


ISPSTP=NCYSTPMNPTPCY+1 

10900 


I5HSTP=NCYSHSMNPTPCY+1 

11000 


RSHTCT=RSHT 

11100 


RSHTT=ASHTT/RSHTCT 

11200 


AIE=0.0 

11300 


WTWHEL=MTWHELm9.0 

11900 


AIWHEL = (WTWHEL/2,0)}<(RWHELM><M2+(RWHELri/2.0))<X2) 

11500 


AIVEH=WTCAR»RWHELMi<X2 

11600 


CALL ICSTUP 

11700 


STOP 

11800 


END 







s 


TABLE XVI. - SAMPLE OUTPUT FOR SUPPLY TRANSIENT (100 POINTS/CYCLE) 


STIRLING ENGINE FOUR CYLINDER, THREE VOLUMES PER CYLINDER, CONTROLS MODEL 


RUN CONDITIONS FOR THIS TRANSIENT 
TWH = 978.0 TWC = 359.0 


TWH = 978.0 

TWC = 

359.0 

CYCLPR = 

5.000 

NONENG 


1 NUMBCY = 

1000 

ALEAK = 0. 

0000 

ARLEAK 

= 0.0000 

PSOURC = 

10.00 

ISUP5T 


501 I5HTST = 

5000001 

TIME 

POWER! 

PFRICT 

PAUX 

PRRF 

PDRAG 

PLOAD 

PENG 

PNET 

ITRAN 

PAVEMP 

POSDEG 

TORQT 

TFRICT 

TAUX 

TRRF 

TDRAG 

TLOAD 

TENG 

TNET 

KWORK 

SPDRPM 

0 .0000 

0.0000 

3.201 

2.108 

2.531 

2.359 

4.890 

5.309 

-10.20 

1 

4 . 93 

270.0 

52.25 

15.29 

10.07 

12.09 

11.26 

23.35 

25.35 

3.543 

1 

2000.00 

0.3000E-01 

10.67 

3.221 

2.108 

2.531 

2.360 

4.891 

5.329 

0.3469 

101 

5.00 

270 . 0 

48.84 

15.38 

10.07 

12.09 

11.27 

23.35 

25.45 

0.4489E-01 

101 

2000.06 

0.6000E-01 

10.67 

3.234 

2.108 

2.531 

2.360 

4.891 

5.342 

0.3330 

201 

5.00 

270.0 

48.85 

15.44 

10.07 

12.09 

11.27 

23.35 

25.51 

-C.9811E-02 

101 

2000.12 

0.8999E-01 

10.67 

3.298 

2.108 

2.531 

2.360 

4.891 

5.407 

0.2684 

301 

5.00 

270 . 0 

48.86 

15.75 

10.07 

12.09 

11.27 

23.35 

25.82 

-0.3075 

101 

2000.17 

0.1200 

10.67 

3.264 

2.108 

2.531 

2.360 

4.891 

5.372 

0.3022 

401 

5.00 

270 . 0 

48.85 

15.59 

10.07 

12.09 

11.27 

23.35 

25.65 

-0.1534 

101 

2000.23 

0.1500 

10.67 

3.215 

2.108 

2.531 

2.360 

4.892 

5.323 

0.3522 

501 

5 .00 

270.0 

48.99 

15.35 

10.07 

12.09 

11.27 

23.35 

25.41 

0.2215 

101 

2000.29 

0.1800 

11.38 

3.245 

2.108 

2.532 

2.361 

4.892 

5.354 

1.028 

601 

5.05 

270.0 

50.14 

15.49 

10.07 

12.09 

11.27 

23.36 

25.56 

1.221 

101 

2000.45 

0.2100 

11.61 

3.328 

2.109 

2.532 

2.362 

4.893 

5.436 

1.177 

701 

5.17 

270 . 0 

51 . 30 

15.89 

10.07 

12.09 

11.27 

23,36 

25.95 

1.994 

101 

2000.64 

0.2399 

11.84 

3.344 

2.109 

2.532 

2.362 

4.894 

5.453 

1.384 

801 

5.29 

270 . 0 

52.41 

15.96 

10.07 

12.09 

11.27 

23.36 

26.03 

3.023 

101 

2000.85 

0.2699 

12.06 

3.343 

2.109 

2.532 

2.363 

4.896 

5.452 

1.599 

901 

5.40 

270 . 0 

53.45 

15.96 

10.07 

12.09 

11.28 

23.36 

26.02 

4.066 

101 

2001.09 

0.2999 

12.27 

3.395 

2.109 

2.533 

2.364 

4.897 

5.505 

1.753 

1001 

5.51 

270 . 0 

54 . 48 

16.20 

10.07 

12.09 

11.28 

23.37 

26.27 

4.850 

101 

2001.35 

0.3299 

12.47 

3.478 

2.110 

2.533 

2.365 

4.898 

5.588 

1.870 

1101 

5.62 

270.0 

55.50 

16.59 

10.07 

12.09 

11.28 

23.37 

26.66 

5.468 

101 

2001.64 

0.3598 

12.67 

3.498 

2.110 

2.533 

2.366 

4.900 

5.608 

2.041 

1201 

5.72 

270.0 

56.43 

16.69 

10.07 

12.09 

11.29 

23.37 

26.76 

6.305 

101 

2001.94 

0.3898 

12.86 

3.478 

2.110 

2.534 

2.367 

4.901 

5.588 

2.246 

1301 

5.81 

270.0 

57.36 

16.59 

10.07 

12.09 

11.29 

23.38 

26.65 

7.325 

101 

2002.28 

0.4198 

13.04 

3.524 

2.111 

2.534 

2.369 

4.903 

5.635 

2.378 

1401 

5.91 

270 . 0 

58.24 

16.80 

10.07 

12.09 

11.29 

23.38 

26.87 

7.992 

101 

2002.63 


VKPH 

60.56 

60.56 

60.56 

60.56 

60.56 

60.57 

60.57 

60.58 

60.58 

60.59 

60.60 

60.61 

60.62 

60.63 


60.64 



0.^^97 

13.21 

5.609 

2.111 

2.535 

270.0 

59.09 

17.21 

10.07 

12.09 

0.^797 

13.38 

3.631 

2.112 

2.535 

270.0 

59.91 

17.31 

10.07 

12.09 

0.5096 

13.59 

3.595 

2.112 

2.536 

270.0 

60.72 

17.13 

10.07 

12.09 

0.5396 

13.70 

3.697 

2.113 

2.536 

270.0 

61.51 

17.38 

10.07 

12.09 

0.5695 

13.85 

3.729 

2.113 

2.537 

270.0 

62.25 

17.79 

10.07 

12.09 

0 . 599 ^ 

19.00 

3.738 

2.113 

2.538 

270.0 

62.99 

17.81 

10.07 

12.09 

0.6293 

19.19 

3.700 

2.119 

2.538 

270.0 

63.66 

17.62 

10.07 

12.09 

0.6592 

19.28 

3.798 

2.115 

2.539 

270.0 

69.32 

17.89 

10.07 

12.09 

0.6891 

19.91 

3.822 

2.115 

2.539 

270.0 

69.96 

18.19 

10.07 

12.09 

0.7190 

19.53 

3.816 

2.116 

2.590 

270.0 

65.57 

18.16 

10.07 

12.09 

0.7989 

19.65 

3.790 

2.116 

2.591 

270.0 

66.17 

18.03 

10.07 

12.09 

0.7788 

19.77 

3.833 

2.117 

2,591 

270.0 

66.79 

18.23 

10.07 

12.09 

0.8087 

19.88 

3.907 

2.117 

2.592 

270.0 

67.29 

18.57 

10.07 

12.09 

0.8385 

19.99 

3.899 

2.118 

2.593 

270.0 

67.83 

18.51 

10.07 

12.09 

0.8689 

15.10 

3.865 

2.119 

2.599 

270.0 

68.32 

18.37 

10.07 

12.09 

0.8982 

15.20 

3.909 

2.119 

2.599 

270.0 

68.85 

18.57 

10.07 

12.09 

0,9281 

15.30 

3.981 

2.120 

2.595 

270.0 

69.32 

18.90 

10.07 

12.09 

0.9579 

15.39 

3.972 

2.121 

2.596 

270.0 

69.78 

18.85 

10.07 

12.09 

0.9877 

15.98 

3.937 

2.121 

2.597 

270.0 

70.22 

18.69 

10.07 

12.09 


2.370 

11.30 

9.905 

23.38 

5.720 

27.28 

2.969 

8.932 

1501 

101 

6.00 

2003.01 

60.65 

2.371 

11.30 

9.907 

23.39 

5.792 

27.37 

2.608 

9.151 

1601 

101 

6.09 

2003.90 

60.66 

2.373 

11.31 

9.909 

23.39 

5.707 

27.20 

2.802 

10.12 

1701 

101 

6.17 

2003.81 

60.67 

2.379 

11.31 

9.911 

23.90 

5.759 

27.99 

2.903 

10.67 

1801 

101 

6.25 

2009.29 

60.69 

2.376 

11.32 

9.913 

23.90 

5.837 

27.81 

2.973 

11.03 

1901 

101 

6.33 

2009.69 

60.70 

2.378 

11.32 

9.915 

23.91 

5.852 

27.87 

3.100 

11.66 

2001 

101 

6.90 

2005.15 

60.71 

2.379 

11.33 

9.917 

23.91 

5.819 

27.69 

3.275 

12.56 

2101 

101 

6.98 

2005.64 

60.73 

2,381 

11.33 

9.920 

23.92 

5.862 

27.91 

3.359 

12.99 

2201 

101 

6.55 

2006.13 

60.79 

2.383 

11.39 

9.922 

23.93 

5.938 

28.26 

3.911 

13.27 

2301 

101 

6.61 

2006.65 

60.76 

2.385 

11.35 

9.925 

23.93 

5.931 

28.22 

3.539 

13.92 

2901 

101 

6.68 

2007.17 

60.77 

2.387 

11.35 

9.927 

23.99 

5.906 

28.10 

3.682 

19.69 

2501 

101 

6.79 

2007.71 

60.79 

2.389 

11.36 

9,930 

23.99 

5.950 

28.30 

3.751 

15.00 

2601 

101 

6.80 

2008.27 

60.81 

2.391 

11.36 

9.933 

23.95 

6.029 

28.69 

3.786 

15.20 

2701 

101 

6.86 

2008.89 

60.82 

2.393 

11.37 

9,936 

23.96 

6.013 

28.58 

5.902 

15.79 

2801 

101 

6.92 

2009.91 

60.89 

2.395 

11.38 

9.939 

23.96 

5.989 

28.93 

9.031 

16.92 

2901 

101 

6.97 

2010.01 

60.86 

2.397 

11.38 

9.992 

23.97 

6.029 

28.69 

9.089 

16.79 

3001 

101 

7.03 

2010.61 

60.88 

2.399 

11.39 

9.999 

23.98 

6,101 

28.97 

9.105 

16.87 

3101 

101 

7.08 

2011.23 

60.90 

2.902 

11.90 

9.997 

23.98 

6.092 

28.92 

9.209 

17.37 

3201 

101 

7.12 

2011.85 

60.92 

2.909 

11.91 

9.951 

23.99 

6.059 

28.75 

9.329 

17.97 

3301 

101 

7.17 

2012.99 

60.99 



TABLE XVI 
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CT\ 




1.017 

270.0 

15.57 

70.69 

3.977 

18.87 

2.122 

10.07 

2.598 

12.09 

2.906 

11.91 

1 . 0^7 

270.0 

15.65 

71.05 

9.099 

19.20 

2.123 

10.07 

2.598 

12.09 

2.908 

11.92 

1.077 

270.0 

15.73 

71.95 

9.036 

19.19 

2.123 

10.07 

2.599 

12.09 

2.911 

11.93 

1.107 

270.0 

15.81 

71.81 

3.995 

18.93 

2.129 

10.07 

2.550 

12.09 

2.913 

11.99 

1.136 

270.0 

15.89 

72.21 

9.035 

19.12 

2.125 

10.07 

2.551 

12.09 

2.916 

11.99 

1.166 

270.0 

15.96 

72.56 

9.105 

19.99 

2.126 

10.07 

2.552 

12.09 

2.918 

11.95 

1.196 

270.0 

16.03 

72.91 

9.115 

19.98 

2.126 

10.07 

2.553 

12.09 

2.921 

11.96 

1.226 

270.0 

16.10 

73.23 

9.053 

19.18 

2.127 

10.07 

2.559 

12.09 

2.923 

11.97 

1.255 

270.0 

16.16 

73.55 

9.089 

19.35 

2.128 

10.07 

2.555 

12.09 

2.926 

11.98 

1.285 

270.0 

16.22 

73.85 

9.162 

19.68 

2.129 

10.07 

2.555 

12.09 

2.928 

11.98 

1.315 

270.0 

16.28 

79.12 

9.165 

19.69 

2.130 

10.07 

2.556 

12.09 

2.931 

11.99 

1 . 39 ^ 

270.0 

16.39 

79.93 

9.100 

19.38 

2.130 

10.07 

2.557 

12.09 

2.939 

11.50 

1.379 

270.0 

16.90 

79.70 

9.133 

19.53 

2.131 

10.07 

2.558 

12.09 

2.936 

11.51 

1.909 

270.0 

16.95 

79.95 

9.212 

19.89 

2.132 

10.07 

2.559 

12.09 

2.939 

11.52 

1.933 

270.0 

16.50 

75.21 

9.213 

19.89 

2.133 

10.07 

2.560 

12.09 

2.992 

11.53 

1.963 

270.0 

16,55 

75.95 

9.192 

19.55 

2.139 

10.07 

2.561 

12.09 

2.999 

11.53 

1.992 

270.0 

16.60 

75.66 

9.161 

19.63 

2.139 

10.07 

2.562 

12.09 

2.997 

11.59 

1.522 

270.0 

16.65 

75.92 

9.299 

20.01 

2.135 

10.07 

2.563 

12.09 

2.950 

11.55 


9.959 

23.50 

6.099 

28.93 

9.367 

18.21 

3901 

101 

7.22 

2013.19 

60.95 

9.957 

23.51 

6.172 

29.27 

9.379 

18.28 

3501 

101 

7.26 

2013.79 

60.97 

9.960 

23.51 

6.160 

29.20 

9.963 

18.73 

3601 

101 

7.30 

2019.95 

60.99 

9.963 

23.52 

6.119 

29.00 

9.578 

19.28 

3701 

101 

7.39 

2015.13 

61.02 

9.967 

23.53 

6.160 

29.18 

9.608 

19.50 

3801 

101 

7.38 

2015.81 

61.09 

9.970 

23.59 

6.231 

29.51 

9.606 

19.51 

3901 

101 

7.92 

2016.50 

61.06 

9.973 

23.55 

6.292 

29.55 

9.660 

19.81 

9001 

101 

7.96 

2017.19 

61.08 

9.977 

23.55 

6.180 

29.25 

9.786 

20.93 

9101 

101 

7.99 

2017.90 

61.10 

9.980 

23.56 

6.217 

29.91 

9.810 

20.57 

9201 

101 

7.53 

2018.61 

61.12 

9.989 

23.57 

6.291 

29.75 

9.799 

20.53 

9301 

101 

7.56 

2019.32 

61.19 

9.987 

23.58 

6.299 

29.76 

9.896 

20.79 

9901 

101 

7,59 

2020.09 

61.16 

9.991 

23.59 

6.230 

29.95 

9.965 

21.90 

9501 

101 

7.62 

2020.78 

61.19 

9.995 

23.59 

6.269 

29.59 

9.983 

21.51 

9601 

101 

7.65 

2021.52 

61.21 

9.998 

23.60 

6.399 

29.96 

9.953 

21.39 

9701 

101 

7.68 

2022.26 

61.23 

5.002 

23.61 

6.396 

29.96 

9.999 

21.69 

9801 

101 

7.70 

2023.00 

61.25 

5.005 

23.62 

6.276 

29.62 

5.115 

22.21 

9901 

101 

7,73 

2023.76 

61.28 

5.009 

23.63 

6.295 

29.70 

5.139 

22.39 

5001 

101 

7.76 

2029.52 

61.30 

5.013 

23.69 

6.379 

30.08 

5.098 

22.20 

5101 

101 

7.78 

2025.28 

61.32 



1.551 

16.69 

9.226 

2.136 

270.0 

76.13 

19.92 

10.07 

1.581 

16 . 7 ^ 

9.179 

2.137 

270.0 

76 . 3 <» 

19.69 

10.07 

1.611 

16.78 

9.212 

2.138 

270.0 

76.59 

19.89 

10.07 

1 . 6^0 

16.82 

9.283 

2.139 

270.0 

76.73 

20.17 

10.07 

1.670 

16.86 

9.287 

2.139 

270.0 

76.92 

20.18 

10.07 

1.699 

16.90 

9.219 

2.190 

270.0 

77.10 

19.82 

10.07 

1.729 

16.93 

9.291 

2.191 

270.0 

77.27 

19.95 

10.07 

1.758 

16.97 

9.322 

2.192 

270.0 

77.99 

20.32 

10.07 

1.788 

17.00 

9.296 

2.193 

270.0 

77.60 

20.19 

10.07 

1.817 

17.03 

9.293 

2.199 

270.0 

77.76 

19.93 

10.07 

1 . 8^7 

17.06 

9.260 

2.195 

270.0 

77.89 

20.00 

10.07 

1.876 

17.09 

9.393 

2.196 

270.0 

78.07 

20.39 

10.07 

1.905 

17.12 

9.397 

2.197 

270.0 

78.21 

20.39 

10.07 

1.935 

17.15 

9.270 

2.197 

270.0 

78.39 

20.03 

10.07 

1 . 96 <i 

17.18 

9.286 

2.198 

270.0 

78.97 

20.09 

10.07 

1 . 99 ^ 

17.20 

9.369 

2.199 

270.0 

78.63 

20.97 

10.07 

2.023 

17.23 

9.398 

2.150 

270.0 

78.76 

20.37 

10.07 

2.052 

17.26 

9.295 

2.151 

270.0 

78.88 

20.11 

10.07 

2.082 

17.28 

9.329 

2.152 

270.0 

79.00 

20.29 

10.07 


2.569 

2.953 

5.017 

6.362 

12.09 

11.56 

23.65 

29.99 

2.565 

2.956 

5.020 

6.316 

12.09 

11.57 

23.66 

29.76 

2.566 

2.958 

5.029 

6.399 

12.09 

11.58 

23.66 

29.91 

2.567 

2.961 

5.028 

6.921 

12.09 

11.59 

23.67 

30.29 

2.568 

2.969 

5.032 

6.927 

12.09 

11.60 

23.68 

30.25 

2.569 

2.967 

5.036 

6.359 

12.09 

11.60 

23.69 

29.89 

2.570 

2.970 

5.090 

6.382 

12.09 

11.61 

23.70 

30.02 

2.571 

2,973 

5.099 

6.969 

12.09 

11.62 

23.71 

30.39 

2.572 

2.976 

5.098 

6.939 

12.09 

11.63 

23.72 

30.26 

2.573 

2.979 

5.051 

6,387 

12.09 

11.69 

23.73 

30.00 

2.579 

2.982 

5.055 

6.909 

12.09 

11.65 

23.79 

30.07 

2.575 

2.989 

5.059 

6.989 

12.09 

11.66 

23.75 

30.96 

2.576 

2.987 

5.063 

6.993 

12.09 

11.67 

23.76 

30,96 

2.577 

2.990 

5.067 

6.918 

12.09 

11.68 

23.76 

30.10 

2.578 

2.993 

5.071 

6.939 

12.09 

11.69 

23.77 

30.16 

2.579 

2.996 

5.075 

6.518 

12.09 

11.70 

23.78 

30.59 

2.580 

2.999 

5.080 

6,998 

12.09 

11.71 

23.79 

30.99 

2.581 

2.502 

5.089 

6.996 

12.09 

11.72 

23.80 

30.18 

2.582 

2.506 

5.088 

6.976 

12.09 

11,73 

23.81 

30.31 


5.155 

5201 

7.80 


22.99 

101 

2026.09 

61.35 

5.291 

5301 

7.83 


22.92 

101 

2026.82 

61.37 

5.295 

5901 

7.85 


22.97 

101 

2027.60 

61.39 

5.209 

5501 

7.87 


22.82 

101 

2028.38 

61.92 

5.238 

5601 

7.89 


22.99 

101 

2029.16 

61.99 

5.399 

5701 

7.91 


23.51 

101 

2029.95 

61.96 

5.397 

5801 

7.93 


23.55 

101 

2030.79 

61.99 

5.296 

5901 

7.95 


23.39 

101 

2031.59 

61,51 

5,351 

6001 

7.96 


23.62 

101 

2032.33 

61.59 

5.930 

6101 

7.98 


29.03 

101 

2033.19 

61.56 

5.938 

6201 

8.00 


29.08 

101 

2033.95 

61.58 

5.380 

6301 

8.01 


23.87 

101 

2039.75 

61.61 

5.901 

6901 

8.03 


23.99 

101 

2035.56 

61.63 

5.500 

6501 

8.09 


29.98 

101 

2036.37 

61.66 

5.506 

6601 

8.06 


29.53 

101 

2037.19 

61.68 

5.999 

6701 

8.07 


29.30 

101 

2038.01 

61.71 

5.986 

6801 

8.08 


29.52 

101 

2038.83 

61.73 

5.560 

6901 

8.10 


29.90 

101 

2039.65 

61.76 


5.551 7001 8.11 

2^.88 101 20 < i 0.^8 


61.78 
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^.703 

18.39 

9.888 

2.293 

270.0 

89.70 

22.02 

10.11 

<♦.731 

18 ,< i 0 

9.779 

2.295 

270.0 

8 < i . 7 <V 

21.52 

10.11 

<♦.760 

18.90 

9.670 

2.296 

270.0 

89.77 

21.02 

10.11 

<♦.788 

18.91 

9.700 

2.297 

270.0 

89.82 

21.15 

10.11 

<^.816 

18.92 

9.787 

2.298 

270.0 

89.86 

21.53 

10.11 


18.92 

9.781 

2.299 

270.0 

89.90 

21.99 

10.11 

<♦.872 

18.93 

9.683 

2.250 

270.0 

89.99 

21.05 

10.11 

<♦.901 

18.93 

9.710 

2.251 

270.0 

89.99 

21.16 

10.11 

9.929 

18.99 

9.800 

2.252 

270.0 

85.03 

21.55 

10.11 

<♦.957 

18.99 

9.795 

2.253 

270.0 

85.07 

21.52 

10.11 

<♦,985 

18.95 

9.695 

2.259 

270.0 

85.11 

21.06 

10.11 

5.013 

18.96 

9.722 

2.255 

270.0 

85.15 

21.18 

10.11 

5 . 0 <il 

18.96 

9.810 

2.256 

270.0 

85.19 

21.56 

10.11 

5.069 

18.97 

9.806 

2.257 

270.0 

85.23 

21.53 

10.11 

5.0 98 

18.97 

9.707 

2.258 

270.0 

85.27 

21.08 

10.12 

10.00 

18.98 

5.188 

2.962 

270.0 

91.08 

21.65 

10.28 

10.03 

18.98 

5.162 

2.963 

270.0 

91.09 

21.53 

10.28 

10.05 

18.98 

5.173 

2.965 

270.0 

91.12 

21,57 

10.28 


2.683 

2.810 

5.993 

7.131 

12.09 

12.66 

29.79 

32.13 

2.689 

2.813 

5.997 

7.023 

12.09 

12.67 

29.75 

31.63 

2.685 

2.817 

5.502 

6.915 

12.09 

12.68 

29.76 

31.13 

2.686 

2.820 

5.507 

6.997 

12.09 

12.69 

29.77 

31.26 

2.687 

2.829 

5.511 

7.035 

12.09 

12.70 

29.79 

31.69 

2.688 

2,828 

5.516 

7.030 

12.09 

12.71 

29.80 

31.60 

2.690 

2.831 

5.521 

6.933 

12.09 

12.72 

29.81 

31.16 

2.691 

2.835 

5.525 

6.961 

12,09 

12.73 

29.82 

31.27 

2.692 

2.838 

5.530 

7.052 

12.09 

12.79 

29.83 

31.67 

2.693 

2.892 

5.535 

7.098 

12.09 

12.75 

29.89 

31.63 

2.699 

2.895 

5.539 

6.999 

12.09 

12.76 

29.85 

31.18 

2.695 

2.899 

5.599 

6.978 

12.09 

12.77 

29.86 

31.29 

2.696 

2.852 

5.599 

7.066 

12.09 

12.78 

29.87 

31.67 

2.697 

2.856 

5.553 

7.063 

12.09 

12.80 

29.88 

31.65 

2.699 

2.860 

5.558 

6.966 

12.09 

12.81 

29.89 

31.20 

2.896 

3.535 

6.932 

7,650 

12.09 

19.75 

26.89 

31.92 

2.897 

3.539 

6.937 

7.625 

12.09 

19.76 

26.85 

31.81 

2.898 

3,593 

6.991 

7.637 

12.09 

19,77 

26.86 

31.85 


5.590 

16101 

8.68 


27.83 

101 

2119.95 

69.19 

5.698 

16201 

8,69 


28.36 

101 

2120.83 

69.22 

5.807 

16301 

8.69 


28.88 

101 

2121.72 

69.29 

5.779 

16901 

8.69 


28.7 9 

101 

2122.62 

CM 

5.691 

16501 

8.70 


28.99 

101 

2123.51 

69.30 

5.697 

16601 

8.70 


28.50 

101 

2129.39 

69.32 

5.799 

16701 

8.70 


28.98 

101 

2125.29 

69.35 

5.767 

16801 

8.71 


28.90 

101 

2126.18 

69.38 

5.677 

16901 

8.71 


28.53 

101 

2127.07 

69.90 

5.683 

17001 

8.71 


28.60 

101 

2127.95 

69.93 

5.789 

17101 

8.72 


29.08 

101 

2128.85 

69.96 

5.756 

17201 

8.72 


29.00 

101 

2129.79 

69.99 

5.667 

17301 

8.72 


28.65 

101 

2130.63 

69.51 

5.671 

17901 

8.73 


28.70 

101 

2131.51 

69.59 

5.769 

17501 

8.73 


29.18 

101 

2132.91 

69.57 

9.703 

35601 

9.06 


32.32 

101 

2288.65 

69.30 

9.729 

35701 

9.06 


32.99 

101 

2289.98 

69.32 


< V .710 35801 9.06 

32.91 101 2290.30 


69.35 



10 . OS 
270.0 

1 S . 9 S 

91.11 

5 . 24 S 

21.88 

2.466 

10.28 

2.899 

12.09 

10.11 

270.0 

18 . 9 S 
91 . OS 

5.261 

21.92 

2.467 

10.28 

2.900 

12.09 

10.13 

270.0 

18 . 9 S 
91.11 

5.183 

21.59 

2.468 

10.28 

2.901 

12.09 

10.16 

270.0 

18.99 

91.09 

5.198 

21.64 

2.469 

10.28 

2.903 

12.09 

10 . IS 
270.0 

18.99 
91 . OS 

5.199 

21.64 

2.470 

10.28 

2.904 

12.09 

10.21 

270.0 

18.99 

91.10 

5.210 

21.68 

2.472 

10.28 

2.905 

12.09 

10.24 

270.0 

18. *9 9 
91 . OS 

5.252 

21.85 

2.473 

10.29 

2.906 

12.09 

10.26 

270.0 

18.99 

91.11 

5.193 

21.59 

2.474 

10.29 

2.907 

12.09 

10.29 

270.0 

18.99 

91.09 

5.207 

21.64 

2.475 

10.29 

2.908 

12.09 

10.31 

270.0 

18.99 

91.09 

5.194 

21.58 

2.476 

10.29 

2.909 

12.09 

10.34 

270.0 

18,99 

91.11 

5.202 

21.61 

2.477 

10.29 

2.910 

12.09 

10.37 

270.0 

19.00 

91.09 

5.255 

21.82 

2.478 

10.29 

2.911 

12.09 

14.94 

270.0 

19.12 

S 7.97 

5.637 

22.06 

2.690 

10.53 

3.088 

12.09 

14.96 

270,0 

19.12 

87.94 

5.650 

22.11 

2.691 

10.53 

3.089 

12.09 

14.99 

270.0 

19.12 
S 7. 91 

5.583 

21.84 

2.693 

10.53 

3.090 

12.09 

15.01 

270.0 

19.12 

S 7.88 

5.597 

21.89 

2.694 

10.53 

3.091 

12.09 

15.04 

270.0 

19.12 
87 . S 5 

5.577 

21.80 

2.695 

10.54 

3.092 

12.09 

15.06 

270.0 

19.12 

87.82 

5.591 

21.85 

2.696 

10.54 

3.093 

12.09 

15.09 

270.0 

19.12 

S 7.79 

5.636 

22.02 

2.697 

10.54 

3.094 

12.09 


3.547 

14.78 

6.446 

26.87 

7.714 

32.16 

4.629 

32.08 

35901 

101 

9.06 

2291.12 

69.37 

3.551 

14.79 

6.451 

26.88 

7.728 

32.20 

4.612 

32.00 

36001 

101 

9.06 

2291.94 

69.40 

3.554 

14.80 

6.456 

26.89 

7.651 

31.87 

4.685 

32.35 

36101 

101 

9.06 

2292.77 

69.42 

3.558 

14.81 

6.461 

26.90 

7.667 

31.93 

4.665 

32.27 

36201 

101 

9.07 

2293.59 

69.45 

3.562 

14.83 

6.466 

26.91 

7.669 

31.92 

4.660 

32.25 

36301 

101 

9.07 

2294.41 

69.47 

3.566 

14.84 

6.471 

26.92 

7.681 

31.96 

4.644 

32.21 

36401 

101 

9.07 

2295.24 

69.50 

3.570 

14.85 

6.475 

26.93 

7.725 

32.13 

4.597 

32.01 

36501 

101 

9.07 

2296.05 

69.52 

3.574 

14.86 

6.480 

26.94 

7.667 

31.88 

4.651 

32.28 

36601 

101 

9.07 

2296.87 

69.55 

3.577 

14.87 

6.485 

26.95 

7.682 

31.93 

4.633 

32.21 

36701 

101 

9.07 

2297.69 

69.57 

3.581 

14.88 

6.490 

26.96 

7.670 

31.87 

4.641 

32.25 

36801 

101 

9.07 

2298.52 

69.60 

3.585 

14.89 

6.495 

26.98 

7.679 

31.90 

4.628 

32.24 

36901 

101 

9.07 

2299.34 

69.62 

3.589 

14.90 

6.500 

26.99 

7.734 

32.11 

4.570 

31.99 

37001 

101 

9.07 

2300.15 

69.65 

4.285 

16.77 

7.373 

28.85 

8.327 

32.59 

3.237 

26.52 

55101 

101 

9.20 

2440.15 

73.88 

4.289 

16.78 

7.378 

28.86 

8.342 

32.64 

3.218 

26.44 

55201 

101 

9.20 

2440.87 

73.91 

4.293 

16.79 

7.382 

28.87 

8.276 

32.37 

3.279 

26.66 

55301 

101 

9.20 

2441.60 

73. 93 

4.296 

16.80 

7.387 

28.88 

8.291 

32.42 

3.259 

26.57 

55401 

101 

9.20 

2442.33 

73.95 

4.300 

16.81 

7.392 

28.89 

8.272 

32.34 

3.274 

26.62 

55501 

101 

9.20 

2443.06 

73.97 

4.304 

16.82 

7.397 

28.90 

8.287 

32.39 

3.254 

26.53 

55601 

101 

9.20 

2443.79 

73.99 

4.308 

16.83 

7.401 

28.91 

8.333 

32.56 

3.203 

26.31 

55701 

101 

9.20 

2444.51 

74.02 



TABLE XVI. - Continued 


15.11 

270.0 


15.14 

270.0 


15.16 

270.0 


15.18 

270.0 


15.21 

270.0 


15.23 

270.0 


15.26 

270.0 


15.28 

270.0 


19.81 

269.9 


19.83 

269.9 


19.86 

269.9 


19.88 

269.9 


19.90 

269.9 


19.93 

269.9 


19.95 

269.9 


19.97 

269.9 


19.99 

269.9 


20.02 

269.9 


19.12 

87.76 

5.595 

21.86 

19.12 

87.73 

5.617 

21.93 

19.12 

87.70 

5.589 

21.82 

19.12 

87.67 

5.653 

22.06 

19.12 

87.64 

5.666 

22.11 

19.12 

87.62 

5.600 

21.85 

19.12 

87.58 

5.614 

21.89 

19.12 

87.55 

5.593 

21.81 

19.13 

82.50 

5.890 

21.95 

19.13 

82.49 

5.937 

22.12 

19.13 

82.47 

5.895 

21.96 

19.13 

82.45 

5.915 

22.03 

19.13 

82.43 

5.887 

21.92 

19.13 

82.42 

5.942 

22.12 

19.13 

82.40 

5.961 

22.19 

19.13 

82.39 

5.900 

21.95 

19.13 

82.37 

5.914 

22.00 

19.13 

82.35 

5.891 

21.91 


2.698 

10.54 

3.094 

12.09 

2.699 

10.54 

3.095 

12.09 

2.700 

10.54 

3.096 

12.09 

2.701 

10.54 

3.097 

12.09 

2.702 

10.54 

3.097 

12.09 

2.703 

10.55 

3.098 

12.09 

2.704 

10.55 

3.099 

12.09 

2.705 

10.55 

3.100 

12.09 

2.900 

10.80 

3.244 

12.09 

2.901 

10.81 

3.244 

12.09 

2.902 

10.81 

3.245 

12.09 

2.902 

10.81 

3.246 

12.09 

2.903 

10.81 

3.246 

12.09 

2.904 

10.81 

3.247 

12.09 

2.905 

10.81 

3.248 

12.09 

2.906 

10.81 

3.248 

12.09 

2.907 

10.81 

3.249 

12,09 

2.908 

10.82 

3.250 

12.09 


4.311 

16.84 

7.406 

28.92 

4.315 

16.85 

7.410 

28.93 

4.318 

16.85 

7.414 

28.94 

4.321 

16.86 

7.418 

28.95 

4.324 

16.87 

7.422 

28.96 

4.327 

16.88 

7.426 

28.97 

4.331 

16.89 

7.430 

28.97 

4.334 

16.90 

7.434 

28.98 

4.966 

18.50 

8.210 

30.59 

4.969 

18.51 

8.213 

30.60 

4.972 

18.52 

8.217 

30.60 

4.975 

18.52 

8.221 

30.61 

4.978 

18.53 

8.224 

30.62 

4.981 

18.54 

8.228 

30.63 

4.984 

18.55 

8.232 

30.63 

4.987 

18.55 

8.236 

30.64 

4.990 

18.56 

8.239 

30.65 

4.993 

18.57 

8.243 

30.66 


8.294 

32.39 

3.239 

26.44 

8.316 

32.47 

3.213 

26.32 

8.289 

32.36 

3.236 

26.40 

8.354 

32.61 

3.167 

26.11 

8.369 

32.66 

3.149 

26.03 

8.303 

32.39 

3.210 

26.26 

8.318 

32.44 

3.192 

26.16 

8.299 

32.36 

3.208 

26.21 

8.790 

32.75 

1.958 

19.16 

8.837 

32.92 

1.907 

18.97 

8.797 

32.77 

1.943 

19.10 

8.817 

32.83 

1.919 

19.01 

8.790 

32.73 

1.942 

19.09 

8.846 

32.93 

1.883 

18.86 

8.867 

33.00 

1.859 

18.77 

8.806 

32.77 

1.916 

18.98 

8.821 

32.82 

1.897 

18.91 

8.799 

32.73 

1.915 

18.97 


55801 

101 

9.20 

2445.14 

55901 

101 

9.20 

2445.75 

56001 

101 

9.20 

2446.36 

56101 

101 

9.20 

2446.96 

56201 

101 

9.20 

2447.57 

56301 

101 

9.20 

2448.18 

56401 

101 

9.20 

2448.80 

56501 

101 

9.20 

2449.41 

76101 

101 

9.27 

2563.11 

76201 

101 

9.27 

2563.64 

76301 

101 

9.27 

2564.17 

76401 

101 

9.27 

2564.69 

76501 

101 

9.27 

2565.22 

76601 

101 

9.27 

2565.74 

76701 

101 

9.27 

2566.26 

76801 

101 

9.27 

2566.78 

76901 

101 

9.27 

2567.31 

77001 

101 

9.27 

2567.83 


74.04 

74.05 
74.07 
74.09 
74.11 
74.13 

74.15 

74.16 

77.61 

77.62 
77.64 

77.66 

77.67 

77.69 

77.70 

77.72 

77.73 


77.75 



20 . 0 <i 

269.9 

19.13 

82.33 

5.903 

21.95 

2.909 

10.82 

3.250 

12.09 

20.06 

269.9 

19.13 

82.32 

5.969 

22.12 

2.910 

10.82 

3.251 

12.09 

20.09 

269.9 

19.13 

82.30 

5.908 

21.96 

2.911 

10.82 

3.252 

12.09 

20.11 

269.9 

19.13 

82.29 

5.927 

22.03 

2.912 

10.82 

3.252 

12.09 

20.13 

269.9 

19.13 

82.27 

5.900 

21.92 

2.913 

10.82 

3.253 

12.09 

20.15 

269.9 

19.13 

82.25 

5.956 

22.12 

2.916 

10.82 

3.256 

12.09 

20.18 

269.9 

19.13 

82.23 

5.976 

22.19 

2.915 

10.83 

3.256 

12.09 

20.20 

269.9 

19.13 

82.21 

5.913 

21.96 

2.916 

10.83 

3.255 

12.09 

20.22 

269.9 

19.13 

82.19 

5.927 

22.00 

2.917 

10.83 

3.255 

12.09 

20.25 

269.9 

19.13 

82.18 

5.906 

21.92 

2.918 

10.83 

3.256 

12.09 

20.27 

269.9 

19.13 

82.15 

5.916 

21.95 

2,918 

10.83 

3.257 

12.09 

20.29 

269.9 

19.13 

82.13 

5.962 - 
22.12 

2.919 

10.83 

3.257 

12.09 

20.32 

269.9 

19.13 

82.11 

5.921 

21.97 

2.920 

10.83 

3.258 

12.09 

20.36 

269.9 

19.13 

82.09 

5.960 

22.03 

2.921 

10.83 

3.259 

12.09 

20.36 

269.9 

19.13 

82.07 

5.913 

21.92 

2.922 

10.86 

3.259 

12.09 

20.38 

269.9 

19.13 

82.05 

5.965 

22.12 

2.923 

10.86 

3.260 

12.09 

20.61 

269.9 

19.13 

82.03 

5.986 

22.19 

2.926 

10.84 

3.261 

12.09 

20.63 

269.9 

19.13 

82.02 

5.926 

21.96 

2.925 

10.86 

3.261 

12.09 

20.65 

269.9 

19.13 

82.00 

5.960 

22.01 

2.926 

10.86 

3.262 

12.09 


6.996 

8.267 

8.812 

1.898 

77101 

9.27 


18.58 

30.66 

32.77 

18.90 

101 

2568.36 

77.77 

5.000 

8.250 

8.859 

1.867 

77201 

9.27 


18.58 

30.67 

32.96 

18.71 

101 

2568.88 

77.78 

5.003 

8.256 

8.819 

1.886 

77301 

9.27 


18.59 

30.68 

32.78 

18.86 

101 

2569.60 

77.80 

5.006 

8.258 

8.839 

1.860 

77601 

9.27 


18.60 

30.69 

32.85 

18.75 

101 

2569.92 

77.81 

5.009 

8.262 

8.813 

1.883 

77501 

9.27 


18.61 

30.69 

32.76 

18.83 

101 

2570.66 

77.83 

5.012 

8.265 

8.867 

1.825 

77601 

9.27 


18.61 

30.70 

32.96 

18.61 

101 

2570.96 

77.86 

5.015 

8.269 

8.888 

1.800 

77701 

9.27 


18.62 

30.71 

33.01 

18.51 

101 

2571.66 

77.86 

5.018 

8.273 

8.828 

1.857 

77801 

9.27 


18.63 

30.72 

32.78 

18.71 

101 

2571.98 

77.88 

5.021 

8.276 

8.866 

1.838 

77901 

9.27 


18.66 

30.72 

32.83 

18.66 

101 

2572.50 

77.89 

5.026 

8.280 

8.822 

1.856 

78001 

9.27 


18.66 

30.73 

32.76 

18.70 

101 

2573.01 

77.91 

5.027 

8.283 

8.836 

1.860 

78101 

9.27 


18.65 

30.76 

32.79 

18.63 

101 

2573.53 

77.92 

5.030 

8.287 

8.881 

1.789 

78201 

9.27 


18.66 

30.75 

32.95 

18,63 

101 

2576.03 

77.96 

5.033 

8.291 

8.862 

1.825 

78301 

9.27 


18.67 

30.75 

32.80 

18.56 

101 

2576.55 

77.95 

5.036 

8.296 

8.861 

1.802 

78601 

9.27 


18.67 

30.76 

32.87 

18.66 

101 

2575.05 

77.97 

5.039 

8.298 

8.835 

1.825 

78501 

9.27 


18.68 

30.77 

32.76 

18.56 

101 

2575.56 

77.98 

5.062 

8.302 

8.888 

1.768 

78601 

9.27 


18.69 

30.78 

32.95 

18.33 

101 

2576.07 

78.00 

5.065 

8.305 

8.910 

1.763 

78701 

9.27 


18.70 

30.78 

33.02 

18.23 

101 

2576.57 

78.01 

5.068 

8.309 

8.851 

1.799 

78801 

9.28 


18.70 

30.79 

32.80 

18.63 

101 

2577.08 

78.03 

5.051 

8.312 

8.866 

1.780 

78901 

9.28 


18.71 

30.80 

32.85 

18.36 

101 

2577.58 

78.05 



20.48 

269.9 

19.13 

81.98 

5.917 

21.92 

2.927 

10.84 

3.263 

12.09 

20.50 

269.9 

19.13 

81.96 

5.928 

21.96 

2.928 

10.84 

3.263 

12.09 

20.52 

269.9 

19.13 

81.94 

5.975 

22.12 

2.929 

10.84 

3.264 

12.09 

20.54 

269.9 

19.13 

81.92 

5.934 

21.97 

2.929 

10.85 

3.264 

12.09 

20.57 

269.9 

19.13 

81.91 

5.953 

22.03 

2.930 

10.85 

3.265 

12.09 

20.59 

269.9 

19.13 

81.90 

5.925 

21.93 

2.931 

10.85 

3.266 

12.09 

20.61 

269.9 

19.13 

81.89 

5.977 

22.11 

2.932 

10.85 

3.266 

12.09 

20.64 

269.9 

19.13 

81.87 

5.998 

22.19 

2.933 

10.85 

3.267 

12.09 

20.66 

269.9 

19.13 

81.86 

5.938 

21.96 

2.934 

10.85 

3.268 

12.09 

20.68 

269.9 

19.13 

81.84 

5.952 

22.01 

2.935 

10.85 

3.268 

12.09 

20.70 

269.9 

19.13 

81.83 

5.930 

21.92 

2.936 

10.85 

3.269 

12.09 

20.73 

269.9 

19.13 

81.81 

5.941 

21.96 

2.937 

10.86 

3.270 

12.09 

20.75 

269.9 

19.13 

81.80 

5.987 

22.13 

2.938 

10.86 

3.270 

12.09 

20.77 

269.9 

19.13 

81.78 

5.947 

21.97 

2.939 

10.86 

3.271 

12.09 


OVER TIME LIMIT 


CPU TIMES 


XVI. - Concluded 


5.053 

8.316 

8.844 

1.799 

79001 

9.28 


18.72 

30.80 

32.76 

18.41 

101 

2578.08 

78.06 

5.056 

8.320 

8.856 

1.783 

79101 

9.28 


18.73 

30.81 

32.80 

18.34 

101 

2578.59 

78.08 

5.059 

8.323 

8.903 

1.732 

79201 

9.28 


18.73 

30.82 

32.97 

18.15 

101 

2579.09 

78.09 

5.062 

8.327 

8.863 

1.768 

79301 

9,28 


18.74 

30.83 

32.82 

18.28 

101 

2579.59 

78.11 

5.065 

8.330 

8.883 

1.745 

79401 

9.28 


18.75 

30.83 

32.88 

18.19 

101 

2580.09 

78.12 

5.068 

8.334 

8.857 

1.768 

79501 

9.28 


18.75 

30.84 

32.78 

18.28 

101 

2580.59 

78.14 

5.071 

8.338 

8.909 

1.712 

79601 

9.28 


18.76 

30.85 

32.96 

18.08 

101 

2581.09 

78.15 

5.074 

8.341 

8.931 

1.687 

79701 

9.28 


18.77 

30.86 

33.04 

17.98 

101 

2581.60 

78.17 

5.077 

8.345 

8.872 

1.742 

79801 

9.28 


18.78 

30.86 

32.82 

18.18 

101 

2582.10 

78.18 

5.080 

8.348 

8.887 

1.723 

79901 

9.28 


18.78 

30.87 

32.87 

18.11 

101 

2582.60 

78.20 

5.083 

8.352 

8.865 

1.742 

80001 

9.28 


18.79 

30.88 

32.78 

18.17 

101 

2583.10 

78.21 

5.086 

8.356 

8.878 

1.726 

80101 

9.28 


18.80 

30.88 

32.82 

18.11 

101 

2583.60 

78.23 

5.089 

8.359 

8.924 

1.676 

80201 

9.28 


18 . 81 

30.89 

32.98 

17.92 

101 

2584.10 

78.24 

5.092 

8.363 

8.885 

1.711 

80301 

9.28 


18 .81 

30.90 

32.83 

18.05 

101 

2584.61 

78.26 


20.00 MINUTES. 



TABLE XVII. - SAMPLE OUTPUT FOP SUPPLY TRANSIENT (25 POINTS/CYCLE) 


STIRLING ENGINE FOUR CYLINDER, THREE VOLUMES PER CYLINDER CONTROLS MODEL 


RUN 

CONDITIONS FOR 

THIS TRANSIENT 



TUJH = 

978.0 

TWC = 

359.0 

CYCLPR = 

5.000 

ALEAK = 

0.0000 

ARLEAK 

= 0.0000 

PSOURC = 

10.00 

TIME 

POWERT 

PFRICT 

PAUX 

PRRF 

PDRAG 

POSDEG 

TORQT 

TFRICT 

TAUX 

TRRF 

TDRAG 

0.0000 

0.0000 

3.201 

2.108 

2.531 

2.359 

270.0 

52.25 

15.29 

10.07 

12.09 

11.26 

0.1200 

10.55 

3.214 

2.108 

2.531 

2.360 

270.2 

48.78 

15.35 

10.07 

12.09 

11.27 

0.2400 

11.71 

3.310 

2.109 

2.532 

2.362 

270.2 

52.68 

15.80 

10.07 

12.09 

11.27 

0.3599 

12.52 

3.464 

2.110 

2.534 

2.366 

270.2 

56.56 

16.52 

10.07 

12.09 

11.29 

0.4797 

13.23 

3.573 

2.112 

2.535 

2.372 

270.2 

59.92 

17.03 

10.07 

12.09 

11.30 

0.5995 

13.34 

3.676 

2.114 

2.538 

2.378 

270.2 

62.85 

17.51 

10.07 

12.09 

11.32 

0.7191 

14.36 

3.783 

2.116 

2.540 

2.385 

270.2 

65.40 

18.00 

10.07 

12.09 

11.35 

0.8386 

14.82 

3.849 

2.118 

2.543 

2.393 

270.2 

67.62 

18.29 

10.07 

12.09 

11.37 

0.9579 

15.21 

3.923 

2.121 

2.546 

2.402 

270.2 

69.55 

18.62 

10.07 

12.09 

11.40 

1.077 

15.54 

3.983 

2.124 

2.549 

2.411 

270.2 

71.22 

18.88 

10.07 

12.09 

11.43 

1.196 

15.84 

4.038 

2.127 

2.553 

2.421 

270.2 

72.68 

19.12 

10.07 

12.09 

11.46 

1.315 

16.09 

4.087 

2.130 

2.557 

2.432 

270.2 

73.93 

19.32 

10.07 

12.09 

11.49 

1.434 

16.30 

4.130 

2.133 

2.560 

2.442 

270.2 

75.04 

19.49 

10.07 

12.09 

11.53 

1.552 

16.49 

4.168 

2.136 

2.564 

2.454 

270.2 

75.97 

19.65 

10.07 

12.09 

11.56 

1.670 

16.65 

4.206 

2.140 

2.568 

2.465 

270.2 

76.82 

19.79 

10.07 

12.09 

11.60 


NONENG 

s 

1 NUMBCY = 

1000 


ISUPST 

= 

126 ISHT5T = 

1250001 


PLOAD 

PENG 

PNET 


ITRAN 

PAVEMP 


TLOAD 

TENG 

TNET 


KUORK 

5PDRPM 

VKPH 

4.890 

5.309 

-10.20 


1 

4 . 93 


23.35 

25.35 

3.543 


1 

2000.00 

60.56 

4.892 

5.322 

-0 .8021E- 

01 

101 

5.00 


23.35 

25,41 

0.1758E- 

01 

26 

2000.25 

60.56 

4.895 

5.419 

0.9541 


201 

5.29 


23.36 

25.87 

3.450 


26 

2000.89 

60.58 

4.900 

5.574 

1.576 


301 

5.72 


23.37 

26.59 

6.593 


26 

2002.01 

60.62 

4.907 

5.685 

2.133 


401 

6.09 


23.39 

27,10 

9.436 


26 

2003.49 

60.66 

4.916 

5.790 

2.601 


501 

6.40 


23.41 

27,58 

11.86 


26 

2005.26 

60.72 

4 . 925 

5.898 

2.986 


601 

6.68 


23.43 

28.06 

13.90 


26 

2007.30 

60.78 

4.936 

5.968 

3.340 


701 

6.91 


23.46 

28.36 

15.80 


26 

2009.56 

60.85 

4.948 

6.044 

3.627 


801 

7.12 


23.49 

28.69 

17.37 


26 

2012.01 

60.92 

4.961 

6.106 

3.874 


901 

7.30 


23.52 

28.95 

18.75 


26 

2014.63 

61.00 

4 . 974 

6.165 

4.082 


1001 

7.45 


23.55 

29.18 

19.95 


26 

2017.38 

61.08 

4.988 

6.216 

4.256 


1101 

7.58 


23.58 

29.39 

20 . 96 


26 

2020.26 

61.17 

5.003 

6.263 

4.403 


1201 

7.70 


23.61 

29.56 

21.87 


26 

2023.23 

61.26 

5.018 

6.305 

4.525 


1301 

7.80 


23.65 

29.72 

22.61 


26 

2026.29 

61.35 

5.033 

6.345 

4.623 


1401 

7.88 


23.68 

29.86 

23.27 


26 

2029.42 

61.45 



TABLE XVII. - Continued 


1.789 

270.2 

16.79 

77.52 

4.257 

20.00 

2.143 

10.07 

2.572 

12.09 

2.477 

11.64 

1.906 

270.2 

16.91 

78.10 

4.261 

19.99 

2.147 

10.07 

2.576 

12.09 

2.488 

11.67 

2.02^ 

270.2 

17.02 

78.67 

4.286 

20.08 

2.150 

10.07 

2.580 

12.09 

2.501 

11.71 

2.142 

270.2 

17.11 

79.20 

4.310 

20.15 

2.154 

10.07 

2.585 

12.09 

2.513 

11.75 

2.259 

270.2 

17.20 

79.66 

4.374 

20.42 

2.158 

10.07 

2.589 

12.09 

2.525 

11.79 

2.376 

270.2 

17.28 

80.14 

4.352 

20.28 

2.162 

10.07 

2.593 

12.09 

2.538 

11.83 

2.493 

270,2 

17.36 

80.54 

4.376 

20.36 

2.165 

10.08 

2.598 

12.09 

2.550 

11.87 

2.610 

270.2 

17.42 

80.94 

4.392 

20.40 

2.169 

10.08 

2.602 

12.09 

2.563 

11.90 

2.727 

270.2 

17.49 

81.29 

4.409 

20.45 

2.173 

10.08 

2.606 

12.09 

2.576 

11.94 

2.843 

270.2 

17.54 

81.63 

4.430 

20.51 

2.177 

10.08 

2.611 

12.09 

2.589 

11.98 

2.959 

270.2 

17.60 

81.92 

4,443 

20,53 

2.181 

10.08 

2.615 

12.09 

2.602 

12.03 

3.075 

270.2 

17.65 

82.21 

4.459 

20.58 

2.185 

10.08 

2.619 

12.09 

2.615 

12.07 

3.191 

270.2 

17.70 

82.51 

4.476 

20.62 

2,189 

10.08 

2.624 

12.09 

2.629 

12.11 

3.307 

270.2 

17.74 

82.83 

4.516 

20.77 

2.193 

10.08 

2.628 

12.09 

2.642 

12.15 

3.422 

270.2 

17.79 

83.11 

4,527 

20.78 

2,197 

10.09 

2.633 

12.09 

2.656 

12.19 

3.537 

270.2 

17.83 

83.37 

4.521 

20.72 

2.201 

10.09 

2.637 

12.09 

2.669 

12.23 

3,652 

270.2 

17,87 

83.67 

4.558 

20.85 

2.205 

10.09 

2.642 

12.09 

2,683 

12.27 


ITERATION FAILURE 13 21 

DEBUG OUTPUT FROM FWS3V2, FLOW CONTINUES 


5.049 

6.400 

4.685 

1501 

7.96 


23.72 

30.07 

23.73 

26 

2032.61 

61.54 

5.065 

6.408 

4.775 

1601 

8.02 


23.76 

30.06 

24.28 

26 

2035.85 

61.64 

5.081 

6.437 

4.832 

1701 

8.08 


23.80 

30.15 

24.73 

26 

2039.13 

61.74 

5.097 

6.464 

4.881 

1801 

8.13 


23.83 

30.23 

25.14 

26 

2042.45 

61.84 

5.114 

6.532 

4.879 

1901 

8.17 


23.87 

30.49 

25.30 

26 

2045.80 

61.94 

5.131 

6.514 

4.957 

2001 

8.21 


23.91 

30.36 

25.87 

26 

2049.18 

62.05 

5.148 

6.541 

4.983 

2101 

8.25 


23.95 

30.43 

26.16 

26 

2052.59 

62.15 

5.165 

6.561 

5.011 

2201 

8.29 


23.99 

30.48 

26.47 

26 

2056.01 

62.25 

5.182 

6.582 

5.032 

2301 

8.32 


24.03 

30.52 

26.74 

26 

2059.47 

62.36 

5.200 

6.607 

5.046 

2401 

8.35 


24.07 

30.59 

26.97 

26 

2062.93 

62.46 

5.217 

6.624 

5.063 

2501 

8.37 


24.11 

30.62 

27.19 

26 

2066.41 

62.57 

5.235 

6.644 

5.075 

2601 

8.40 


24.15 

30.66 

27.40 

26 

2069.90 

62.67 

5.253 

6.665 

5.082 

2701 

8.43 


24.19 

30.70 

27.62 

26 

2073.41 

62.78 

5.271 

6.709 

5.063 

2801 

8.45 


24.23 

30.85 

27.74 

26 

2076.93 

62.89 

5.288 

6.724 

5.072 

2901 

8,47 


24.28 

30.87 

27.96 

26 

2080.45 

62.99 

5.306 

6.722 

5.093 

3001 

8.49 


24.32 

30.81 

28.24 

26 

2083,99 

63.10 

5.325 

6.764 

5,069 

3101 

8.51 


24.36 

30.94 

28.37 

26 

2087.53 

63.21 



DEBUG PRINTOUT FROM FWS3V2 NUMBER I 

-1.8^68522 0.67807871 1.7777653 

8.<i^73810 8.^34<»292 8.391955^ 

10.520517 10.571616 10.618118 

8.3913279 8.38^6169 8.<il37^<*0 

6.79<t8027 6.7637110 6.7301712 

933.05029 670.09106 ^07.1318^ 

933.05029 670.09106 ^07.1318^ 

218.9^711 3.63002^0 2^.396973 

0.00000000 0.00000000 0.00000000 

0.00000000 0.00000000 0.00000000 

6.0000000 4.0000000 0.00000000 

DEBUG PRINTOUT FROM FW53V2 NUMBER 2 

1 0.34085358E-03-0.15308842E-05 

2 0.35189348E-03 0 . 1 96686 78E-05 

3 0.42028003E-03-0.36629757E-07 

4 0.26040501E-03 0 . 90656 331 E-0 5 

5 0.44105924E-03-0.21163014E-05 

6 0.41166646E-03-0.36774009E-05 

7 0.15074215E-03-0.11467416E-03 

8 0.34981524E-03 0 .52733361E-04 

9 0.61267894E-03-0.80659611E-05 

10 0.23896278E-03-0.59822432E-05 

11 0.28218934E-03 0 . 13787368E-04 

12 0.58980775E-03-0.39224651E-05 

13 2.7060757 0 . 36422080E-06 

14 218.94711 0.56013025E-07 

3.767 17.90 4.549 2,209 

270.2 83.94 20.78 10.09 


-1.0093956 
0.68167634E-02 
0.16364049E-01 
0.35321091E-02- 
■0.26894122E-01- 
933.05029 
933.05029 
90.201843 
0.00000000 
0 .00000000 
0.00000000 


0.22354629E-01 

0.17652512E-01 

0.15330065E-01 

0.24474800E-01 

670.09106 

670.09106 

34.976685 

0.00000000 

0.00000000 

0.00000000 


0.18807023E 
0, 15282376E 
0.38049284E 
0.10422764E 
0. 17917191E 
0.47990270E 
0 .81206042E 
0.11413758E 
0.63257085E 
0.99459339E 
0.88281422E 
0.54192342E 
0.82929730E 
15.669580 
2.646 
12.09 


-05-0 
-05-0 
05 0 
05 0 
-05-0 
-05-0 
06-0 
05 0 
-05-0 
-06-0 
-06-0 
-05 0 
02 

2.697 

12.31 


.68167634E-02 

.15537865E-01 

.22354629E-01 

.26894122E-01 

.24193227E-02 

.24474800E-01 

.35321091E-02 

.18862173E-01 

.15330065E-01 

.16364049E-01 

.12884624E-02 

.17652512E-01 

218.94711 

3.6300240 

5.343 
24.40 


155.32368 

95.280518 

69.150574 

135.37712 

407.13184 

407.13184 

9.6353846 

0.00000000 

0.00000000 


84.119629 

65.120804 

122.31100 

147.19943 


0.00000000 

0.00000000 


0 .00000000 
0 .00000000 


0.34576468E* 
0.36992598E' 
0.39727031E- 
0.22953982E' 
0 .44039311E* 
0.44319732E 
0.15737723E 
0.32948749E 
0.62636589E' 
0.25784690E' 
0.28549996E 
0.56761317E 
2.4547462 
218.94292 
6.758 
30.87 


03-0.17407069E-02- 
■03-0.15878677E-01- 
■03 0.17734930E-01 
•03 0.26885588E-01 
•03 0.35782866E-02 
-03-0.30463874E-01- 
03-0 .80595501E-02 
•03 0.16584896E-01 
03-0.85253455E-02 
•03-0.16540229E-01- 
03-0.44727325E-02- 
•03 0.21012962E-01 
218.94292 
3.7022800 
5.091 3201 

28,67 26 


0.42787343E-02 

0.15708271E-01 

0.20044778E-01 

0.26889853E-01 

0.57948194E-03 

■0.27469337E-01 

■0.57958290E-02 

0.17723534E-01 

'0.11927705E-01 

-0.16452137E-01 

■0.28805975E-02 

0.19332737E-01 

218.94495 

3.6661520 

8.53 

2091.08 63.31 


ITERATION FAILURE 13 21 

DEBUG OUTPUT FROM FWS3V2. FLOW CONTINUES 
DEBUG PRINTOUT FROM FWS3V2 NUMBER 1 


-1.8488550 0.67807108 1.7777700 

8.4564724 8.4434786 8.4009447 

10.532256 10.583419 10.630010 

8.4006414 8.3939114 8.4231005 

6.8018503 6.7706938 6.7370911 

933.04053 670.08960 407.13867 

933.04053 670.08960 407.13867 

219.13310 3.6251764 24.417908 

0.00000000 0.00000000 0.00000000 

0.00000000 0.00000000 0.00000000 

6.0000000 4.0000000 0.00000000 

DEBUG PRINTOUT FROM FWS3V2 NUMBER 2 

1 0.34122425E-03-0.30216620E-06 
0.35227183E-03 0 . 13753242E-05 
0. 42 0724 04E- 03-0 .256 138 04 E-0 6 
0.26069884E-03 0 . 10765216E-05 
0. 4415526 1 E-03-0. 220817 95E-05 
0.41212025E-03 0 . 75434764E-06 
0. 15091080E-03-0.10253876E-03 
0.35020383E-03 0. 44601446E-04 
0 .61334926 E-03-0. 597 14303E-05 
0.2 3 92 1287 E-0 3-0 . 90203696E-06 
0.28248131E-03 0 . 13715817E-05 

12 0.59040473E-03-0.51222283E-07 

13 2.7060795 0 . 33993888E-06 

14 219.13310 0.75821980E-08 

ITERATION FAILURE 13 21 

DEBUG OUTPUT FROM FW53V2» FLOW CONTINUES 
DEBUG PRINTOUT FROM FWS3V2 NUMBER 1 


2 

3 

4 

5 

6 

7 

8 
9 

10 

11 


-1.8487787 

8.4610968 

10.537929 

8.4049740 

6.8053713 

933.03564 

933.03564 

219.22597 

0.00000000 

0.00000000 

6.0000000 


0.67828947 
8.4481249 
10.589179 
8.3982439 
6.7741909 
670.08887 
670.08887 
3.6283350 
0.00000000 
0 . 00000000 
4.0000000 


1.7776604 

8.4055243 

10.635817 

8.4274426 

6.7405472 

407.14209 

407.14209 

24.428360 

0.00000000 

0.00000000 

0.00000000 


-1.0093880 
0 .68388470E- 
0 . 16398180E- 
0.35421478E- 
-0 .26927751E- 
933.04053 
933.04053 
90.293518 
0 .00000000 
0 .00000000 
0 .00000000 


0.22386253E-01 
0.17685641E-01 
0.15362691E-01 
0.24521481E-01 
670 . 08960 
670.08960 
34.929977 
0.00000000 
0.00000000 
0.00000000 


02 

01 

02 - 

01 - 


0 . 18827 64 lE-05-0. 6838847 OE-02 
0 . 152 98 938 E- 05-0 . 155474 06 E-0 1 
0 . 38089847E-05 0 . 22386 253E-0 1 
0 . 10434396E-05 0 . 26 927 751 E-01 
0 . 17 937 127 E-05-0 .240626 93E-02 
0 .48 04 3 148 E-0 5-0 .2452 148 IE- 01 
0 .81312339E-06-0.35421478E-02 
0.11427246E-05 0 . 18 904S39E-01 
0 .633301 99E-05-0 . 153626 91 E-01 
0 .99563658E-06-0.16398180E-01 
0 .88372758E-06-0 . 1287 460 3E-02 
0.54247266E-05 0 . 17685641E-01 
0.79968683E-02 219.13310 

15.682890 3.6251764 


-1.0095730 
0 .68273023E- 
0.16410727E- 
0.35421478E- 
-0 .26973929E- 
933.03564 
933.03564 
90.333954 
0 .00000000 
0.00000000 
0.00000000 


0.22421386E-01 

0.17707225E-01 

0.15367709E-01 

0.24546076E-01 

670.08887 

670.08887 

34.960403 

0.00000000 

0.00000000 

0.00000000 


02 

01 

02 - 

01 - 


155.32376 
95.280701 
69.150452 
135.37694 
407.13867 
407.13867 
9.6353846 
0 .00000000 
0 .00000000 


84.119797 

65.120743 

122,31084 

147.19954 


0 .00000000 
0.00000000 


0 .00000000 
0.00000000 


0 . 346 14024 E-0 3-0. 17 3367 98E-02-0 .428626 31 E-02 
0 .370 32505E-03-0. 15932381 E-0 1-0. 1573 98 92E-01 
0 . 39768941E-03 0 . 1778066 9E-01 0 . 20 08346 1 E-0 1 
0 .22979786E-03 0 . 26956864E-01 0 . 26 94230 5E-0 1 
0 .44088671E-03 0 . 35657 398 E-02 0 . 57 97 3526 E-0 3 
0. 4436867 3E-03-0 .3 05226 04 E-0 1-0 ,27522 042E-01 
0 . 1 57 5536 3E-0 3-0 .806 958 97 E-0 2-0 .580 5868 7 E-0 2 
0 . 32985443E-03 0 . 166 05478 E-0 1 0 . 17 7551 58E- 0 1 
0 .627 0506 4 E-0 3-0. 85358880 E-02-0 . 1 1 94 92 90 E-0 1 
0.2581 1652 E-03-0. 1656 5826 E-0 1-0 . 164820 03E-01 
0 .2857 9566 E-0 3-0. 4491 3068E-02-0 .288 938 36 £-02 
0 .56818686E-03 0 . 210571 33E-01 0 . 1 937 1 387 E-0 1 
2.4547501 219.12891 219.13098 

219.12891 3.6948719 3.6600237 


155.32195 
95.275513 
69.153061 
135.38133 
407 .14209 
407 . 14209 
9.6353846 
0 .00000000 
0.00000000 


84.115616 

65.122467 

122.31578 

147.19705 


O.QOQOOOOO 

0.00000000 


0.00000000 

0.00000000 



-4 

S» 


TABLE XVII. - Continued. 


DEBUG PRINTOUT FROM FU53V2 NUMBER 2 

1 0,3<*1<»0865E-03 0.223^8377E- 

2 0.352^6601E-03-0.15709^18E- 

3 0.92092893E-03-0.15726909E- 
^ 0.26082693E-03-0.50635179E- 

5 0 . 99179335E-03 0,33051956E- 

6 0.91235285E-03-0.31190107E- 

7 0.15099508E’’03-0.10750998E- 

8 0 . 35038997E-03 0.93569258E- 

9 0.61368523E-03-0.59550219E- 

10 0.2393957^E-03-0.30926918E- 

11 0.28262753E-03 0,90099668E- 

12 0.59069279E-03-0.27902615E- 


3.882 

270.2 


13 
19 
17.99 
89.18 


2.7059689 
219.22597 

9.568 
20.83 


0.33995917E' 

0.180978Q8E- 

2.219 

10.09 


05 0.18837527E-05-0.68273023E-02 

06 0.15307351E-05-0.15599089E-01 

05 0.38108865E-05 0 . 22921 386 E-0 1 

06 0.10939926E-05 0 . 26 97 3 92 9E-0 1 
05 0.17996668E-05-0.29278536E-02 
05 0.98070529E-05-0.29596076E-01 
03 0.81359502E-06-0.35921978E-02 
09 0.11933212E-05 0 . 18909857E-01 
05 0.63369050E-05-0.15367709E-01 
05 0 .99618319E*“06-0 . 16910727E-01 
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12.09 
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12.36 


5.361 
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0.22990997E- 
0.99112327E- 
0.99393935E- 
0.15769908E- 
0.33002986E- 
0.62739039E- 
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0 .28599967E- 
0.56896230E- 
2.9596399 
219.22177 
6.781 
30.92 


03-0.17977339E-02- 
03-0.15913819E-01- 
03 0.17775830E-01 
03 0.26959879E-01 
03 0.35993680E-02 
03-0.30559292E-01- 
03-0.80866553E-02- 
03 0.16633585E-01 
03-0.85969298E-02- 
03~0.16585902E-01- 
03-0.99857822E-02- 
03 0.21071189E-01 
219.22177 
3.6988306 
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28.82 
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6.983 
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9501 
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12 
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12 
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26 
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270 

.2 

89 . 
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12 

.09 
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36 
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.2 
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12 
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.2 
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21.58 
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12 
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59 

18 
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1.2 

89 

.74 

21.90 

10.51 

12 

.09 


3.040 

5.794 

7.148 

4.746 

5601 

8.83 


13.34 

25.43 

31.37 

30.16 

26 

2176.38 

65.90 

3.055 

5.814 

7.164 

4.724 

5701 

8.84 


13.38 

25.47 

31.39 

30.21 

26 

2179.90 

66.00 

3.070 

5.833 

7.179 

4.700 

5801 

8.85 


13.43 

25.51 

31.40 

30.26 

26 

2183.41 

66.11 

3.537 

6.434 

7.665 

3.840 

8901 

9.01 


14.76 

26.84 

31.98 

29.94 

26 

2288.98 

69.31 

3.552 

6.453 

7.680 

3.809 

9001 

9.01 


14.80 

26.88 

32.00 

29.89 

26 

2292.27 

69.41 

3.567 

6.472 

7.643 

3.831 

9101 

9.01 


14.84 

26.93 

31.80 

30.08 

26 

2295.54 

69.51 

3.583 

6.492 

7.682 

3.776 

9201 

9.02 


14.88 

26.97 

31.91 

29.94 

26 

2298.81 

69.60 

3.598 

6.511 

7.695 

3.747 

9301 

9.02 


14.92 

27.01 

31.92 

29.92 

26 

2302.08 

69.70 

3.613 

6.531 

7.708 

3.718 

9401 

9.02 


14.97 

27.05 

31.93 

29.88 

26 

2305.33 

69.80 

3.629 

6.550 

7.721 

3.689 

9501 

9.03 


15.01 

27.10 

31.94 

29.86 

26 

2308.58 

69.90 

3.644 

6.569 

7.709 

3.685 

9601 

9.03 


15.05 

27.14 

31.85 

29.96 

26 

2311.81 

70.00 

3.659 

6.589 

7,747 

3.630 

9701 

9.03 


15.09 

27.18 

31.96 

29.83 

26 

2315.04 

70.10 

3.674 

6.608 

7.782 

3.579 

9801 

9.04 


15.14 

27.22 

32.06 

29.72 

26 

2318.26 

70.19 

3.690 

6.627 

7.753 

3.592 

9901 

9.04 


15.18 

27.26 

31.90 

29.88 

26 

2321.47 

70.29 

3.705 

6.647 

7.762 

3.566 

10001 

9.04 


15.22 

27.31 

31.89 

29.91 

26 

2324.68 

70.39 

3.720 

6.666 

7.777 

3.534 

10101 

9.05 


15.26 

27.35 

31.91 

29.88 

26 

2327.88 

70,48 

3.736 

6.686 

7.790 

3.504 

10201 

9.05 


15.30 

27.39 

31.92 

29.88 

26 

2331.06 

70.58 

3.751 

6.705 

7.803 

3.474 

10301 

9.05 


15.34 

27.43 

31.92 

29.88 

26 

2334.24 

70.68 

4.225 

7.299 

8.240 

2.485 

13401 

9.12 


16.61 

28.70 

32,40 

28.64 

26 

2428.74 

73.54 
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1^.78 

270.2 

18.79 

89.70 

5.583 

21.90 

2.681 

10.52 

3.081 

12.09 

^.256 

16.69 

7.337 

28.78 

8.264 
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26 

9.13 

2434,57 
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32.60 

25.28 

26 

2563.62 

77.62 

8.231 

8.761 

1.053 

18601 

9.21 


30.63 

32.61 

25.21 

26 

2566.09 

77.70 

8.248 

8.771 

1.025 

18701 

9.21 


30.67 

32.61 

25.14 

26 

2568.54 

77.77 

8.265 

8.782 

0.9974 

18801 

9.21 


30.70 

32.62 

25.07 

26 

2570.97 

77.85 

8.283 

8.792 

0.9692 

18901 

9.21 


30.74 

32.63 

25.01 

26 

2573.41 

77.92 

8.300 

8.803 

0,9412 

19001 

9.21 


30.77 

32.64 

24.95 

26 

2575.83 

77.99 

8.317 

8.813 

0.9135 

19101 

9.21 


30.81 

32.65 

24.88 

26 

2578.26 

78.07 

8.335 

8.831 

0.8780 

19201 

9.21 


30.84 

32.68 

24.83 

26 

2580.67 

78.14 

8.352 

8.849 

0.8431 

19301 

9.22 


30.88 

32.72 

24.75 

26 

2583.06 

78.21 

9.124 

9,307 

- 0,3788 

23901 

9.28 


32.42 

33.08 

21.92 

26 

2687.30 

81.37 

9.140 

9.319 

- 0.4072 

24001 

9.28 


32.46 

33.09 

21.88 

26 

2689.43 

81.43 

9.157 

9.329 

- 0.4332 

24101 

9.28 


32.49 

33.10 

21.86 

26 

2691.57 

81.50 

9.173 

9.332 

- 0.4534 

24201 

9.28 


32.52 

33.09 

21.85 

26 

2693.70 

81.56 

9.189 

9.346 

- 0.4822 

24301 

9,28 


32.55 

33.11 

21.82 

26 

2695.83 

81.63 

9.206 

9.351 

- 0.5041 

24401 

9.28 


32.59 

33.10 

21.81 

26 

2697.95 

81.69 

9.222 

9.378 

- 0.5472 

24501 

9.28 


32.62 

33.17 

21.72 

26 

2700.06 

81.75 

9.238 

9,414 

- 0.5995 

24601 

9.29 


32.65 

33.27 

21.59 

26 

2702.17 

81.82 





o\ 


TABLE XVII. - Concluded 


25.06 

270.2 

18.80 

87.53 

6.211 

21.93 

3.169 

11.19 

3.^22 

12.09 

5.832 

20.60 

9.255 

32.68 

9.380 

33.13 

- 0.5825 

21.73 

29701 

26 

9.29 

2709.28 

81.88 

25.15 

270.2 

18.80 

87.58 

6.218 

21.99 

3.179 

11.20 

3.425 

12.09 

5.896 

20.63 

9.271 

32.71 

9.392 

33.19 

- 0.6115 

21.72 

29801 

26 

9.29 

2706.38 

81.95 

25 . 2 <t 

270.2 

18.80 

87.59 

6.228 

21.96 

3.178 

11.21 

3.428 

12.09 

5.860 

20.66 

9.287 

32.75 

9.906 

33.17 

- 0.6920 

21.68 

29901 

26 

9.29 

2708.99 

82.01 

25.33 

270.2 

18.80 

87.62 

6.252 

22.03 

3.182 

11.21 

3.430 

12.09 

5.873 

20.69 

9.309 

32.78 

9.939 

33.29 

- 0.6865 

21.61 

25001 

26 

9.29 

2710.58 

82.07 


terminated: STOP 

CPU TIME= 13.^8 MINUTES. 


TABLE XVIII. - INPUT ROUTINE (MAINSE) FOR SHORT-CIRCUIT TRANSIENT 


100 C INPUT DATA FOR THE FOUR WORKING SPACE THREE VOLUME 

200 C STIRLING ENGINE MODEL. 

300 C 

^00 COMMON/PST/ AD, AR » RODL *VOE,VOC, VR>RER. RRC» STROKE, ALPHA tPSIMD 

500 COMMON/CARLOD/GTRAN,GR,AIE,AIWHEL,AIVEH,RWHEEL,AINERT 

600 COMMON/CONST/R,G,PIE,DEGR 

700 C0MM0N/MATXIN/VDELTA,FRAC,T0L1,T0L2.T0LSS.N,NTMAX,MPAS, 

800 1 TOLPCG,REF 

900 COMMON/SWITCH/ISS.ICALC, MATRIX, NONENG,IPRTOP,IBRYTH, 

1000 1 IHPCHV, NOBUG 

1100 COMMON /SHRTCT/ WOl IN,WD2IN,WD3IN,WDAIN,WD10UT, WD20UT,WD30UT,WD40tJ 

1200 IT,WDSHT.MSTMXD,WSTMHD,PSHTMX,P5HTMN,WSMX1,WSMN1,VSH,RSHT,RSHTT 

1300 COMMON/RUNCON/ALEAK, ARLEAK,PSOURC,CYCLPR,SPDRPM,POSDEO, 

1600 1 mstotl,twhin,twcin,cyclpm,spdmax 

1500 COMMON/CYDATA/NPTPCY, ITRMAX,NUMBCY 

1600 C0MM0N/TRANDS/I5UPST,I5HTST,TIME.KW0RK,ISPSTP,ISHSTP 

1700 COMMON/OUTPT/P05DEC,WTOT,ITRAN,UT1,WT2,WT3,WT6,WD1,WD2,WD3,WD6, 

1800 1WD5.WD6,WD7,WD8,TORQ1,TORQ2,TORQ3,TORQ^,POWER1,POWER2,POWER3, 

1900 2POWER6,POUERT,TORQT,TLOAD,TENG,PLOAD,PEN6,PNET,PAVEMP,CCCHK 

2000 C 

2100 C 

2200 C ENGINE GEOMETRY 

2300 AD=23.7613 

2600 AR=1.1290 

2500 RQDL^IO.O 

2600 VR=115.339 

2700 V0E=63.87 

2800 VOC=61.70 

2900 STROKE=^.00 

3000 ALPHA=90.0 

3100 C 

3200 C HEATER AND COOLER WALL TEMPERATURES 

3300 TWHIN=705.0 

3A00 TWCIH=86.0 

3500 C 

3600 C FLOW RESISTANCES BETWEEN VOLUMES 

3700 RER=1.9 

3800 RRC=1.9 

3900 C 

<♦000 C CONSTANTS 

6100 DEGR=57.296 

6200 PIE=3.1A16 

6300 R=9125.6 

6600 G=10017.0 

9500 C 



^600 

C 

MATRIX CONVERGENCE INPUT 

^700 


VDELTA=.01 

4800 


FRAC=1.0 

4900 


TOL1=.001 

5000 


TOL2=.01 

5100 


TOL5S=.0001 

5200 


N = 14 

5300 


NTMAX=16 

5400 


MPAS=20 

5500 


T0LPCG=.5 

5600 

c 


5700 

c 

SWITCHES 

5800 


ISS = 0 

5900 


ICALC^l 

6000 


MATRIX=1 

6100 


NONEHG=l 

6200 


IPRTOP=100 

6300 


IBRYTH=0 

6400 


IHPCNV=0 

6500 


NOBUG=l 

6600 

c 


6700 

c 

RUN CONDITIONS 

6800 


ALEAK=0.0 

6900 


ARLEAK=0.0 

7000 


PSOURC=10.0 

7100 


CYCLPM=15.0 

7200 


SPOMAX=4000.0 

7300 


CYCLPR=15.0 

7400 


SPDRPM=4000.0 

7500 


POSDEO=270.0 

7600 


WSTOTL=. 007786 

7700 


VSH=32.78 

7800 


RSHT=38.0 

7900 


ASHTT=.20 

8000 

c 


8100 

c 

LOAD CHARACTERISTICS 

8200 


GTRAN=1.0/2.53 

8300 


GRS2.84 

8400 


WTENG-0.0 

8500 


UTWHEL=29.48 

8600 


WTCAR=1420.0 

8700 


RWHEEL=30.48 

8800 

c 


8900 

c 

CYCLE DATA 

9000 


NPTPCY=25 

9100 


NUMBCY=1000 

9200 

c 


9300 

c 

TRANSIENT DESIRED 

9400 


NCYSHT=5 

9500 


NCYSUP=50000 

9600 


NCYSTP=50000 

9700 


NCYSHS=60000 

9800 

c 


9900 

c 

CALCULATED INPUT 

10000 


PSHTMX=CYCLPR^1.0 

10100 


PSHTMN=CYCLPR-1.0 

10200 


IP (NONENG ,EQ. 0) N^12 

10300 


REF=(TOL1*TOL2)/2.0 

10400 


RWHELM^RWHEEL/IOO.O 

10500 


ITRMAX=HPTPCYXNUMBCY+1 

10600 


ISUPST-KCYSUPMNPTPCY+l 

10700 


ISHTST=NCYSHT«NPTPCY+1 

10800 


ISPSTP=NCYSTP»NPTPCY+1 

10900 


ISHSTP=NCYSHS*<NPTPCY+1 

11000 


RSHTT=ASHTT/RSHT 

11100 


AIE=0.0 

11200 


WTWHEL=WTWHEL*<4.0 

11300 


AIWHEL = (WTWHEL/2.0)»(RUHELMJ<»2+(RWHELM/2.0)NM2) 

11400 


AIVEH=WTCARmRUHELMNM2 

11500 


CALL ICSTUP 

11600 


STOP 

11700 


END 



TABLE XIX 


UPDATE STRATEGIES FOR 


JACOBIAN MATRIX 


Broyden 

algorithm 

TOLPCG 

Strategy 

used 



No 

1.0 - 0.5 

Use current matrix 


<0.5 

Generate new matrix 

Yes 

1.0 - 0.7 

Use current inverse 


0.7 - 0.0 

Update inverse with algorithm 


<0.0 

Generate new matrix 
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